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Filing date 
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Unit 

Examiner 



09/826,118 
01/09/2001 

Urbain A. von der Embse 
2124 

Chat C. Do 



02/14/2006 



Response to Office Action 02/07/2006 

I have corrected the amendment to the specification in the 
substitute specification by underlining the new paragraph in page 
9 lines 12-35, to be compliant to the Office Action. 



ENCLOSURES 

• PTO/SB/21 Transmittal form 

• A copy of the Notice of Non-Compliant Amendment (37 CFR1.121) 
Office Action mailed to me on 02/07/2006 

• Marked up version of the substitute specification with the 
underlining of the new paragraph on page 9 lines 12-35 to be 
compliant with the "Notice of Non-Compliant Amendment" 



Thanks for all of your help and guidance. 
Sincerely, 

310.641.0488 
Urbain A. von der Embse 
7323 W. 85 th St. 
Westchester, CA 90045-2444 



Contact No. 
Address 



Signature 
Name 



Urbain A. von der Embse 
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Serial No. 



09/826,118 



06/29/2006 



Filing date 



01/09/2001 



Name 



Urbain A. von der Embse 



Unit 



2193 



Examiner 



Chat C. Do 



Response to Office Action 04/10/2006 



PTO-1449 information disclosure is submitted. The abstract 
is amended to respond to your guidelines and is limited to 144 
words. Claims 1-6 are amended to respond to your objections and 
to follow your recommendations. The substitute specification is 
the original specification with editing corrections. Amended 
drawings are the original drawings with editing corrections that 
follow your guidelines. 



• PTO/SB/21 Transmittal form 

• Copy of the "Office Action" mailed to me on 04/10/2006 

• PTO/SB/08a (replaces PTO-1449) information disclosure 

• PTO/SB/08b (replaces PTO-1449) information disclosure 

• Marked up version of the amended abstract of the disclosure 

• Clean version of how the abstract of the disclosure will read 

• Complete list of claims with their status identifiers and 
current amendments with strikethrough and underline markups. 

• Clean version of how the claims will read 

• Marked up version of the substitute specification 

• Clean version of how the specification will read 

• Amended drawings with proper identification in the top margins 



ENCLOSURES 
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REMARKS 

1 . Office Action 

The amendments being submitted are in reference to the 
amendment filed on 02/21/2006. 

2. Office Action 

Claims 1-6 are amended to address your objections and 
recommendations . 

3. Information Disclosure Statement 

References are listed on form PTO/SB/08a for patents and on 
form PTO/SB/08b for technical articles which forms are 
replacements for form PTO-144 9. 

4 . Drawings 

Figure 1-2 are designated as "Prior Art" in their legends. 
Corrected drawings are labeled as "Deleted Sheet" in their top 
margins and replacement sheets are labeled as "Replacement 
Sheet" in their top margins. 

5 . Drawings 

Corrected original sheets are labeled as "Deleted Sheet" 
and corrected sheets are labeled as "Replacement Sheet" in the 
top margins. Drawings 1-8 are the originals. Shading is 
removed from drawings 1-5 and drawing 2 has a correction in 
element 8. Original drawings as part of the specification have 
their cross-out markings removed. 

6 . Drawings 

New drawings 5B-5L have been deleted since they list the 
Matlab software used to calculate performance in the 
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specification and are objected to by the examiner as not being 
supported by the original specification. 

Remaining drawings 1-8 are the original drawings with 
editing corrections and corrections recommended in your "Office 
Action" and are consecutively numbered from 1-8. No drawings 
have been deleted. 

7 . Specification 

The abstract of the disclosure has been amended to follow 
your guidelines. As it now reads the abstract describes the 
enclosure sufficiently to assist the reader in deciding if there 
is a need to consult the full patent in that it identifies the 
combined requirements for a least squares LS solution for 
Wavelets and finite impulse response FIR filters that include 
requirements not currently available in FIR design algorithms, it 
points out that the harmonics of the Wavelets are the design 
coordinates being optimized and these are mapped into the FIR 
time response for implementation, and it points out that these 
harmonic design coordinates provide a single design to generate 
Wavelets at arbitrary parameters which capability is not 
currently available with any design algorithm. 

8 . Specification 

The amended abstract has been limited to 144 words. 

9. Claim Objections 

The referenced "for" and "means for generating complex 
wavelet waveform and filters" have been removed in the amended 
claim 1 which lists the steps required to implement the 
iterative eigenvalue least squares LS algorithm derived in the 
specification. 
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10. Claim Objections 

Objections to claims 3-4 are believed to have been 
corrected in the amended claims 1-6. Claim 1 lists the detailed 
steps in the method claimed for generating the digital mother 
Wavelet at baseband using an iterative eigenvalue least squares 
LS algorithm. Claims 3-4 list the steps and design for 
generating multi-resolution Wavelets from the mother Wavelets in 
claims 1-2. 

11. Claim Rejections - 35 USC § 112 

Requirements for one or more of the claims "particularly 
pointing out and distinctly claiming the subject matter which the 
appliant regards as his invention" are believed to be met by the 
amended claims 1-6. Claims 1-2 list the steps in the method for 
calculating the mother Wavelet time response \|/(n) using the 
design harmonic coordinates and claims 3-5 list the steps 

and designs for the mother Wavelet i|/ k (k) to the generate multi- 
resolution Wavelet time response for scale and translation 
parameters. Claim 6 has been amended to list the properties of 
the Wavelets in claims 1-5. 

12. Claim Rejections - 35 USC § 112 

Objections to claim 1-6 are believed to have been 
corrected in the amended claims which particularly point out and 
distinctly claim the subject matter which is regarded as the 
invention. 

"Scale" is defined in the clean version of the amended 
claim 3 in lines 9-10 in page 5. "Scale" parameters are the 
dilation p, symbol interval M, and symbol length L and wherein p 
can be either positive or negative. "Scale" parameter p in 
current Wavelets is restricted to positive values. 
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Translation parameters q, k are defined in claim 3 in lines 
11-13 in page 5. 

"Fourier harmonics \|/ k (k)" are defined in the clean version 
of the amended claim 1 in lines 4-15 in page 2 in the inverse 
discrete Fourier transform format which derives the time response 
vj/(n) as a function of the frequency response harmonics v|/ k (k). 
Designation of vj/ k ( k) as "design harmonics" is in the clean 
version of the amended claim 3 in line 17 in page 5. 

13. Claim Rejections - 35 USC § 112 

"Means" has been eliminated in the amended claims. 

14. Claim Rejections - 35 USC § 112 

Amended claims 1-6 have been rewritten to eliminate the 
objection that the claims are "narrative in form and replete 
with indefinite and functional or operational language". 

Amended claims 1-2 clearly and positively define the steps 
comprising the method for designing the LS algorithms using 
design coordinates v|/ k (k) to find the optimum y(n) for 
implementation and using one sentence formats. Likewise claims 3- 
5 clearly and positively define the steps and design parameters 
to enable \|/ k (k) in claim 1-2 to generate the multi-resolution 
Wavelet examples. Claim 6 clearly and positively summarizes the 
properties of the new Wavelets which are the subject matter of 
this invention. 
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15,16 Claim Rejections - 35 USC § 101 

Amended claims 1-2 are believed to be in compliance under 
35 USC § 101 since they define in clear detail the steps in the 
new method for designing the mother Wavelet time response \\f (n) 
which is the format required for direct implementation in 
communications and radar applications . 

Amended claims 3-5 are believed to be in compliance under 
35 USC § 101 since they define in clear detail the steps and 
design parameters for generating the multi-resolution Wavelet 
time response for scales p,M,L and sampling rate 1/T from \\i (n) 
using the design harmonics i|/ k (k)for the mother Wavelet and which 
time response is the format required for direct implementation in 
communications and radar applications and which mplementations 
will add the translations q,k. 

Amended claim 6 is in compliance under 35 USC § 101 since 
it lists the properties of the Wavelets for implementation in 
communications and radar transmitters and receivers. 

Demonstrated performance for linear communications, BEM, 
SAR is summarized in the specification drawings in FIG. 4,7,8 
respectively. 

17,18 Claim Rejections - 35 USC § 102 

Amended claims 1-6 are believed to be in compliance under 
35 USC § 102 since the amended claims 1-6 clearly define the 
subject matter of this invention that was not patented or 
described in a printed publication prior to this invention 
application and was not anticipated by prior art summarized in 
pages 2-8 in the specification. 
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In claims 1-2 prior art Wavelet and FIR filter designs 
select time domain coordinates \|/(n) to optimize performance for 
this baseband waveform whereas this invention replaces \j/(n) with 
i|/ k (k), combines Wavelet and FIR requirements, and finds the 
optimum i|/ k (k) to generate the optimum \\f (n) at baseband and for 
all scales p,M,L and all translations q,k. 

Claims 3-5 demonstrate how \j/ k ( k) generates the optimal 
baseband time response i|/(n p ) wherein n p is the sample index for 
scale p, when scale p down-samples or equivalently dilates the 
digital sampling, when scale p dilate M,L at a constant sample 
rate, when scale p up-samples, and by obvious implication for 
any set of scales p,M,L and design parameter 1/T. Current 
Wavelets can only be dilated or dowm-sampled using scaling 
equations or equivalently the trellis construction in FIG. 2. 
Current FIR designs support up-sampling followed by a bandwidth 
limiting filter. 

Application examples in FIG. 6,7,8 demonstrate improved 
performance offered by combining Wavelet and FIR requirements. 

18. Claim Rejections - 35 USC § 102 

Amended claims 1-6 are believed to clearly demonstrate they 
have not been anticipated by prior art in compliance with 35 USC 
§ 102. 

51 of the "Office Action" correctly observes that prior art 
on pages 2-8 provides methods for designing waveforms and filters 
and the Wavelet art in lines 15-35 in page 5 defines the multi- 
resoluion Wavelet as a function of the mother Wavelet at 
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baseband. This definition enables the mother Wavelet to generate 
multi-resolution waveforms at lower bandwidths to tile the 
time-frequency t-f space in Figure 1 using the procedure 
described in lines 6-16 in page 7. There is no mention of 
"using a subset of the Fourier harmonics as the desicn 
coordinates (harmonics) specifying the waveform design" and no 
mention of "design harmonics" in prior Wavelet art and in prior 
waveform and filter art. In line 15 in page 2 to line 4 in page 
4 the waveform and filter design techniques in categories CI to 
C7 including the Wavelet in category C5 are described. There is 
no prior art suggesting the concept of a "means for using the 
design harmonics and frequency translation property to generate 
wavelet waveforms at multiple scales and frequencies". 
"Frequency translation" is a common property of all waveform and 
filter design except possibly category C7. 

52 of the "Office Action" correctly observes that prior art 
in pages 2-8 discloses in lines 15-37 in page 5 the Wavelet 
equation, properties, and its use for generating multi-resolution 
Wavelets at scale p with time translation q which can be used to 
tile the t-f space in Figure 1 as described in lines 6-16 page 7. 
Scale p is a dilation parameter which means the mother Wavelet at 
scale p=0 can generate multi-resolution replicas at scales 
p=l,2,.... corresponding to reductions in sample rate rate 
bandwidth by the factor 2 A p and a stretching by the factor 2 A p to 
tile the t-f space in Figue 1. Another method to tile the t-f 
space uses the iterated filter bank construcion in line 18 page 8 
to line 10 page 9 to generate the multi-resolution Wavelets at 
scales "p=l,2..." with a reduced capability to design the multi- 
resolution Wavelets to be replicas of the mother Wavelet due to 
the multi-stage filtering impacts. Wavelets are real so the 
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multi-resolution Wavelets are real. There is no concept nor use 
of "Fourier harmonics" as design harmonics since all 
implementations are in the time domain as convolutions and sub- 
sampling operations. Frequency response is a performnce metric 
of this design. 

Current Wavelet multi-resolution Wavelet generation is 
restricted to p=l,2,... with corresponding subsampling reducing the 
Nyquist bandwidth 1/T to l/T2 A p. The new Wavelets support the 
design of multi-resolution waveforms with this dilation as well 
as with dilation of bandwidth with constant 1/T, negative values 
for scale p corresponding to up-sampling, and complete 
flexibility in the control of the scale parameters p,M,L. and 
the design sampling rate 1/T. Design coordinates \\f k (k) for the 
mother Wavelet generate all of these multi-resolution waveforms 
as demonstrated in claims 3-5. 

<I2 of the "Office Action" correctly observes that prior art 
in line 15 page 3 to line 4 page 4 includes designing waveforms 
and filters using analytical equations such as the square-root 
raised-cosine definition of a FIR waveform time response that 
allows it to generate multi-resolution waveforms, iterative 
design techniques such as the Remez Exchange and eigenvalue LS 
techniques for designing filters with passband and stopband 
requirements, and polyphase filters which have multi-resolution 
capabilities. There is no concept nor use of "Fourier harmonics" 
as design harmonics since all implementations are in the time 
domain. Frequency response is a performnce metric of these 
designs . 
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19. Response to Amendment 

The Matlab code has been deleted from the drawings and its 
description has been deleted from the specification. Amended 
substitute specification is the original specification with some 
minor editing corrections. Amended drawings are the original 
drawings with some minor edition corrections including your 
requested insertion of "Prior Art" in the legend of drawings 1-2. 
Abstract has been amended to follow your guidelines and consists 
of 144 words. Claims 1-2 have been amended to clearly define the 
method for designing the mother Wavelets described in the 
specification. Claims 3-5 illustrate the derivation of multi- 
resolution Wavelets from the mother Wavelets in claims 1-2 using 
the equations defined in the specification. Claim 6 summarizes 
the claimed properties of the new Wavelets disclosed in claims 1- 
5. 

Objection 1 is believed to be satisfied in that the claims 
only refer to matter in the specification. 

Objection 2 referencing lines 13-35 page 9, Figure 5, and 
Matlab code pages 37-44 is believed to be satisfied in that the 
Matlab code 5.0 and all reference to it has been deleted both 
from the amended drawings and the substitute specification. 

20. Response to Arguments 

The amended claims 1-6 have been rewritten to address your 
objections. Your detailed objections and recommendations were 
most helpful. 
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21. Response to Arguments 

I am extremely grateful for pointing out all of my mistakes 
and lack of understanding of the USPTO rules. 

Thanks ever for all of your help and guidance. 
Sincerely, 

Name Urbain A. von der Embse 

Contact No. 310.641.0488 

Address Urbain A. von der Embse 

7323 W. 85 th St. 

Westchester, CA 90045-2444 
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- 7/ie MAILING DATE of this communication appears on the cover sheet with the correspondence address - 

The amendment document filed on 05 July 2006 is considered non-compliant because it has failed to meet the 
requirements of 37 CFR 1.121 or1.4. In order for the amendment document to be compliant, correction of the following 
item(s) is required.. 

THE FOLLOWING MARKED (X) ITEM(S) CAUSE THE AMENDMENT DOCUMENT TO BE NON-COMPLIANT: 
1. Amendments to the specification: 

□ A. Amended paragraph(s) do not include markings. 

□ B. New paragraph(s) should not be underlined. 
O C. Other See Continuation Sheet 

I3 2. Abstract: 

□ A. Not presented on a separate sheet. 37 CFR 1.72. 
K B. Other See Continuation Sheet 



"New Sheet," or 



3. Amendments to the drawings: 

□ A. The drawings are not properly identified in the top margin as "Replacement Sheet, 

"Annotated Sheer as required by 37 CFR 1.121(d). . 

□ B. The practice of submitting proposed drawing correction has been eliminated. Replacement drawings 

showing amended figures, without markings, in compliance with 37 CFR 1.84 are required. 
G>3 C. Other See Continuation Sheet 



□ 



4. Amendments to the claims: 

□ A. A complete listing of all of the claims is not present. 

□ B The listing of claims does not include the text of all pending claims (including withdrawn claims) 

□ C Each claim has not been provided with the proper status identifier, and as such, the individual status 

of each claim cannot be identified. Note: the status of every claim must be indicated after its claim 
number by using one of the following status identifiers: (Original), (Currently amended), (Canceled), 
(Previously presented), (New), (Not entered), (Withdrawn) and (Withdrawn-currently amended). 

□ D. The claims of this amendment paper have not been presented in ascending numerical order. 
^ E. Other: See Continuation Sheet 

5. Other (e.g., the amendment is unsigned or not signed in accordance with 37 CFR 1.4): 



For further explanation of the amendment format required by 37 CFR 1.121, see MPEP § 714. 
TIME PERIODS FOR FILING A REPLY TO THIS NOTICE: 

1 . Applicant is given no new time period if the non-compliant amendment is an after-final amendment or an amendment 
filed after allowance. If applicant wishes to resubmit the non-compliant after-final amendment with corrections, the 
entire corrected amendment must be resubmitted. 

2 Applicant is given one month, or thirty (30) days, whichever is longer, from the mail date of this notice to supply the 
correction, if the non-compliant amendment is one of the following: a preliminary amendment, a non-final amendment 
(including a submission for a request for continued examination (RCE) under 37 CFR 1 .1 14), a supplemental 
amendment filed within a suspension period under 37 CFR 1.103(a) or (c) t and an amendment filed in response to a 
Quayle action. If any of above boxes 1. to 4. are checked, the correction required is only the corrected section of the 
non-compliant amendment in compliance with 37 CFR 1.121. 

Extensions of time are available under 37 CFR 1.136(a) only if the non-compliant amendment is a non-final 
amendment or an amendment filed in response to a Quayle action. 

Failure to timely respond to this notice will result in: 

Abandonment of the application if the non-compliant amendment is a non-final amendment or an amendment 
filed in response to a Quayle action; or 

Non-entry of the amendment if the non-compliant amendment is a preliminary amendment or supplemental 
amendment. 
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APPLICATION NO. 09/826, 118 

TITLE OF INVENTION: ^few -Wavelet Multi-Resolution Waveforms 
INVENTOR: Urbain Alfred A . von der Embse 



BACKGROUND OF THE INVENTION 



I. Field of the Invention 

TECHNICAL FIELD 

The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data 
communications with data rates up to multiple Tl (1.544 Mbps) and 
higher (>100 Mbps) , and to optical CDMA. Applications are 
mobile, point-to-point and satellite communication networks. More 
specifically the present invention relates to a new and novel 
mcano — for a new approach to the design of waveforms and filters 
using mathematical formulations which generalize the Wavelet 
concept to communications and radar. 
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II. Description of the Related Art 



Multi-resolution waveforms used for signaling and/or 
5 filters and which are addressed in this invention are defined to 
be waveforms of finite extent in time and frequency, with scale 
and shift properties of multi-resolution over the time-frequency 
(t-f) space. These waveforms can also be referred to as multi- 
scale waveforms, and include multi-rate filtering and the 
10 Wavelet as special cases. The emphasis will be on digital 
design and applications with the understanding that these multi- 
resolution waveforms are equally applicable to analog design and 
applications . 

15 Background art consists of the collection of waveform and 

filtering design techniques which can be grouped into six broad 
categories. These categories are: CI) least squares (LS) design 
algorithms for filters and waveforms that design to 
specifications on their frequency response, C2) analytic 

20 filters and waveforms which are specified by a few free design 
parameters that can be sub-categorized into current applications 
and primarily theorectical studies, C3) combinations of CI and 
C2 for greater flexibility in meeting communications and radar 
performance goals, C4) special design techniques to control the 

25 noise levels from intersymbol interference (ISI) and adjacent 
channel interference (ACI) in the presence of timing offsets for 
multiple channel applications, C5) Wavelet filter design using 
scaling functions (iterated filter banks) as the set of design 
coordinates or basis functions, C6) filter and waveform design 

30 techniques for non-linear channels and in particular for 
operation in the non-linear and saturation regions of a high 
power amplifier (HPA) such as a traveling wave tube (TWT) or a 
solid state amplifier, and C7) LS dynamic filters derived from 
discrete filtering and tracking algorithms that include adaptive 

35 equalization for communications, adaptive antenna filters, 
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Wiener filters, Kalman filters, and stochastic optimization 
filters . 

Category CI common examples of LS digital filter design are 
5 the eigenvalue algorithm in "A New Approach to Least-Squares FIR 
Filter Design and Applications Including Nyquist Filters" [7] and 
the Remez-Exchange algorithm in "A Computer Program for Designing 

Optimum FIR Linear Phase Filters' . - E-&-H The eigenvalue algorithm 

is a direct LS minimization and the Remez-Exchange can be 

10 reformulated as an equivalent LS gradient problem through proper 
choice of the cost function. Both LS algorithms use the FIR 
(finite impulse response) digital samples as the set of design 
coordinates. LS design metrics are the error residuals in 

meeting their passband and stopband ideal performance as shown 

15 in FIG. 3. Category C2 common examples of analytical waveforms 
and filters for system applications are the analog Chebyshev, 
Elliptic, Butterworth, and the digital raised-cosine, and square- 
root raised-cosine. For theorectical studies, common examples 
are polyphase multirate filters, quadrature mirror filters 

20 (QMF) , and perfect reconstruction filters. Although these 
theorectical studies have yet to yield realizable useful filters 
for system applications, their importance for this invention lies 
in their identification and application of ideal performance 
metrics for filter designs. Category C3 common digital example 

25 is to start with the derivation of a Remez-Exchange finite 
impulse rooponoo — (-FIR-)- filter and thsen up-sample and filter with 
another bandwidth limiting filter. This results in an FIR over 
the desired frequency band that is larger than available with the 
Remez-Exchange algorithm and with sidelobes that now drop off 

30 with frequency compared to the flat sidelobes of the original 
Remez-Exchange FIR. A category C4 common example is to select 
the free parameters of the category C3 filter in order to 
minimize the signal to noise power ratio of the data symbol (SNR) 
losses from the ISI, ACI, and the non-ideal demodulation. A 

35 second common example is start with a truncated pulse whose 
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length is short enough compared to the symbol repetition 
interval, to accommodate the timing offsets without significant 
impact on the ISI and ACI SNR losses. This shortened pulse can 
then be shaped in the frequency domain. 

5 

Category C5 Wavelet filter design techniques discussed in 
the next section will serve as a useful reference in the 
disclosure of this invention. Category C6 common example is the 
Gaussian minimum shift keying (GMSK) waveform. This is a 

10 constant amplitude phase encoded Gaussian waveform which has no 
sidelobe re-growth through a non-linear or saturating HPA. 
Category C7 common examples are the adaptive equalization 
filter for communication channels, the adaptive antenna filter, 
and the Kalman filter for applications including target tracking 

15 and prediction as well as for equalization and adaptive antennas. 

Minimizing excess bandwith in the waveform and filter 
design is a key goal in the application of the C1,...,C6 design 
techniques for communications and radar. Excess bandwith is 

20 identified as the symbol a in the bandwidth- time product 
BT s =l+a where the two-sided available frequency band is B and the 
symbol repetition interval is T s . Current performance capability 
is represented by the use of the square-root raised cosine (sq-rt 
rc) waveform with a= 0.22 to 0.4 as shown in FIG. 46. The goal 

25 is to design a waveform with <x=0 within the performance 
constraints of ISI, ACI, passband, sideband, and passband 
ripple. This goal of eliminating the excess bandwidth corresponds 
to the symbol rate equal to the available frequency band 1/T S = 
B/(l+a) = B for a=0 . This symbol rate 1/T S =B is well known to be 

30 the maximum possible rate for which orthogonality between 
symbols is maintained. A fundamental performance characteristic 
of our new waveform designs is the ability to eliminate excess 
bandwidth for many applications. 
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Scope of this invention will include all of the waveform 
and filter categories with the exception of category C4 special 
design techniques and category C7 LS dynamic filters. Emphasis 
will be on the category C5 Wavelets to establish the background 
5 art since Wavelets are multi-resolution waveforms that can 
eliminate the excess bandwidth and have known design algorithms 
for FIR waveforms and filters. However, they do not have a design 
mechanism that allows direct control of the ISI, ACI, passband, 
sideband, and passband ripple. Category C2 theorectical studies 

10 also eliminate the excess bandwidth. However, they are not 

multi-resolution waveforms and do not have realizable FIR design 
algorithms. So the emphasis in background art will be on 
Wavelets whose relevant properties we briefly review. 

15 Wavelet background art relevant to this invention consists 

of the discrete Waveform equations and basic properties, 
application of Wavelets to cover a discrete digital time- 
frequency (t-f) signal space, and the design of Wavelets using 
the iterated filter construction. Wavelets are waveforms of 

20 finite extent in time (t) and frequency (f) over the t-f space, 
with multi-resolution, scaling, and translation properties. 
Wavelets over the analog and digital t-f spaces respectively are 
defined by equations (1) and (2) as per Daubechies' s "Ten 
Lectures on Wavelets", Philadelphia : SIAM, 1992 

25 rof crcncc — [1] . 

Continuous wavelet 

= H -,/ V(— ) (i) 

a 

Discrete wavelet 

30 ¥ («) = | a |- ,/ V(— ) (2) 

a 

where the two index parameters "a,b" are the Wavelet dilation and 
translation respectively or equivalently are the scale and shift. 
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The y is the "mother" wavelet and is a real and symmetric 
localized function in the t-f space used to generate the doubly 

i r I/2 

indexed Wavelet vj/ ab . The scale factor " has been chosen to 

keep the norm of the Wavelet invariant under the parameter change 
5 "a,b". Norm is the square root of the energy of the Wavelet 
response. The Wavelets \|; a ,b and \|/ are localized functions in the 
t-f space which means that both their time and frequency lengths 
are bounded. The discrete Wavelet has the time "t" replaced by 
the equivalent digital sample number "n" assuming the waveform is 
10 uniformly sampled at "T" second intervals. 

Wavelets in digital t-f space have an orthogonal basis that 
is obtained by restricting the choice of the parameters "a,b" to 
the values a=2' p , b=qM2 p where *p,q" are the new scale and 
15 translation parameters and "M" is the spacing or repetition 
interval T S =MT of the Wavelets (which from a communications 
viewpoint are symbols) at the same scale "p" . Wavelets at "p,q" 
are related to the mother Wavelet by the equation -ft}- 

20 

^(n) =r^K2- p n-qM) (3) 



25 where the mother Wavelet is a real and even function of the 
sample coordinates. The orthonormality property means that these 
Wavelets satisfy the orthogonality equation with a correlation 
value equal to "1". 

^\|/ pq \|/^ m = 1 iff both p = k and q = m 

n 

30 =0 otherwise (4) 
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Wavelet representation of a digital t-f space starts with 
selecting an N sample time window of a uniform stream of 
digital samples at the rate of 1/T Hz (1/second) equivalent to a 
"T" second sampling interval. The N point or sample t-f space 
5 in FIG. 1 illustrates a Wavelet representation or "tiling" with 
Wavelets that are designed analytically or by an iterated filter 
construction. 

The t-f space in FIG. 1 is partitioned or covered or tiled 
10 by a set of Wavelet subspaces §V p ,p = 0,l,...,m -l) 1 where N=2 m . Each 
Wavelet subspace W p 1 at scale "p" consists of the set of Wavelet 
time translations {q = 0,1,..., N/2 P+1 -1} over this subspace. These 
Wavelet subspaces are mutually orthogonal and the Wavelets within 
each subspace are mutually orthogonal with respect to the time 
15 translates. This N-point t-f space extends over the time 
interval from 0 to (N-l)T 2 whcro — ¥ — is — tke — digital — sampling 
interval, — and over the frequency interval from 0 to (N-l) 3 in 
units of the normalized frequency f NT 4. 

20 The iterated filter bank in FIG. 2 is used to generate the 

Wavelets which cover the t-f space in FIG. 1. Each filter stage 
5 consists of a high pass filter (HPF) and a low pass filter 
(LPF) . Output 6 of the LPF is subsampled by 2 which is 
equivalent to decimation by 2. This t-f space space is an N- 

25 dimensional complex vector metric space V 7. At stage m in the 
iterated filter bank, the remaining t-f space V m _ x 8 is 

partitioned into V m+1 9 and the Wavelet subspace W m+1 10 . 

Scaling functions and Wavelets at each stage of this filter 
30 bank satisfy the following equations 
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^(n)=2- ,/2 Zh q( p(2n-q) Vq 

Kn)=2 ,/2 Zg q ^(2n-q) Vq 
q 



where (p is the scaling function, vy is the Wavelet, HPF P 

coefficients are { h q , Vq } , LPF P coefficients are {g^, Vq}, and 

the equations apply to the stages 0, 1 , ...,m-l . Identifying the 
scale parameter and using the previous Wavelet formulations 
enable these equations to be rewritten for stages p=0, l,...,m-l 
as 



10 < (6) 



For our application the HPF P and LPF P are quadrature mirror 
filters (QMF) with perfect reconstruction- This means they cover 
the subspace V p with flat responses over the subband frequency 
15 including the edges of the frequency subband, and the HPF P 
coefficients are the frequency translated coefficients for the 

LPF P : {g q =(-l) q h q / v ql- 

Wavelet design using iterated filter bank starts with the 
20 selection of the scaling functions. Starting with a primative 
scaling function such as the one proposed by Daubechies — [1] , 
_^one can use the iterated filter construction given by equations 
(5) and (6) to derive successive approximations to a desired 
scaling function which has properties that have been designed 

25 into it by the selection of the filter coefficients {g 9 ,Vq} at 

each level of iteration. The Wavelets can be derived from these 
scaling functions using the iterated filter construction or 
scaling equations (5) and (6) . 
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Another use of the iterated filter construction is to 
design the scaling functions as Wavelets thereupon ending up with 
a larger set of Wavelets for multi-resolution analysis and 
synthesis as illustrated by Coifman's Wavelets in "Wavelet 
5 analysis and signal processing". — Design oxamplco abound in 
the open literature and the recent publications — — H-H — [5] , 
-f-6-J — arc representative illustrations of the divcroc applications 
of Wavelets. 

10 

SUMMARY OF INVENTION 

SUMMARY OF THE INVENTION 

15 This invention is a new and novel design of Wavelet mult -i— 

re solution waveforms and filters in the frequency domain with — a 
pr operty which provides a single waveform design for all of t4*e 
wa veforms at multiple resolutions and scales. Wavelets are- 
d esigned to moot the specific application requirements , are 

20 complex, arc defined by frequency harmonics, include a frequenc y 
tr anslation parameter in addition to the scale and tranolat -iof* 
pa rameters in current Wavelets, and have the property that ^the 
-se t of design coefficients in the frequency domain rent ed** 
jrft variant for all possible multi - resolution scales, and ti ffte— 

25 frequency translations. Frequency design harmonics and the 
f requency translation capability enable the waveform at multiple 
r esolution scales to be derived from the single waveform design 
by scaling the dilation, translation, and frequency translation 
p arameters using the set of frequency design harmonics. Wavelet 

30 waveform design is illustrated by Matlab 5.0 code to design —a 
j ^ncar filter waveform using an iterative least - squa res- 
ei genvalue approach to minimize the non linear least squares c ee£ 
f unction, and to scale this waveform design for multi - rcsolut -ief* 
ap plication specified by the dilation, time translation, -andt 
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•f- rcqucncy tranolation parameters . Additional rcoulto arc given 
le r a conotant amplitude minimum 1 ohift keying bandwidth - cf f icicnt 
modulation waveform and for a synthetic aperture radar waveform. 

5 This invention 4rs — a — sew — approach — to — for the design of 

multi-resolution Wavelet waveforms that improves their 

performance for engineering and scientific applications, and in 
particular for applications to communications and radar. Current 
practice is to 1) design traditional multi-resolution waveforms 

10 in the time domain using metrics which specify frequency 
performance and without consideration of Wavelet properties 7 — and 
or to 2) design Wavelet multi-resolution waveforms using the 
time-domain iterated multi-resolution filtering approach with a 
set of scaling functions used to perform the filtering and 

15 without direct considerations of the frequency performance. 
These two separate design approaches yield fundamentally 
different waveforms. Our ncw This invention provides a means to 
combine these two approaches to generate a new multi-resolution 
Wavelet w aveform with the best properties of the traditional 

20 multi-resolution waveforms and the Wavelet multi-resolution 
waveform. 

Gw^ — This invention introduces ^he — ftew— innovations for the 
design of our new — Wavelet waveform that 1) include the ISI 

25 (intersymbol interference) , ACI (adjacent channel interference) , 
and QMF (quadrature mirror filter) requirements into the design 
algorithm along with non-linear modifications to the traditional 
passband and stopband frequency requirements which are used in 
the Remez-Exchange and eigenvalue LS algorithms— ift — rcf crcncco 

30 [ 8 ] , — — respectively , 2) select the frequency harmonics as the 
design coordinates, 3) use an appropriate subset of the 
available Fourier domain frequency harmonics as the design 
coordinates with the property that this subset of coordinates is 
a basis for the multi-resolution Wavelet waveform design, and 4) 
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design the multi-resolution Wavelet waveform for no excess 
bandwidth a=0. 

These innovations support the development of our new 
5 waveforms 1) that are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar over their t-f space, 2) that can be 
designed to be an orthonormal basis or orthonormal set of 
coordinates over the available frequency bandwidth with the 

10 implicit property that the excess bandwidth vanishes a=0, to 
within the accuracies allowed by communication and radar design 
implementations, and 3) that have multi-scale properties which 
allow a single dc waveform design to be used to uniquely define 
the complete set of now waveforms over the t-f space for multi- 

15 scale applications Specific design algorithm examples developed 
in this invention disclosure are the LS design algorithms that 
use eigenvalue and gradient search techniques in the frequency 
domain to find the best waveform design which minimizes the 
corresponding LS error residual cost function that is a weighted 

20 linear sum of the residual error metrics. 

Wavelet Waveform — waveform performance calculated in this 
invension disclosure from the application of these LS design 
algorithms, demonstrates the capability of this invention to 

25 provide a means to design waveforms that are improved over 
current practice. There are other algorithmic design algorithms 
which can be realized by modifications to this invention. An 
example given in this invention disclosure is the modification of 
these LS algorithms for application to constant amplitude 

30 bandwidth efficient (BEM) waveforms for communications. 
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BRIEF DESCRIPTION OF THE DRAWINGS AND 



THE PERFORMANCE DATA 

BRIEF DESCRIPTION OF DRAWING S AND PERFORMANCE DATA 

5 The above-mentioned and other features, objects, design 

algorithms, and performance advantages of the present invention 
will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 
10 elements, and in which: 

FIG. 1 is a Wavelet a«— N-point t-f space extending over the 
time interval (0, (N-l)T] and the frequency interval (0, (N- 
1)/NT], which is tiled or covered by a set of orthonormal 
15 Wavelets at the scales p=0, 1, ...,m-l . 

FIG. 2 is a Wavelet iterated filter bank used to generate 
the set of Wavelets which tile or cover the t-f space in FIG. 1. 

20 FIG. 3 illustrates the power spectral density (PSD) of our 

new Wavelet waveform at dc, and the set of stopband and passband 
design requirements . 

FIG. — 3 io the power spectral denoity — (PSD) — template for the 
25 the otopband and paosband used to construct the LS error metrics, 



FIG. 4 is the flow diagram of the LS metrics and cost 
functions and the final cost function used to find the optimal LS 
30 solution for our the Wavelet waveform. 

FIG. 5 is flow diagram for the LS recursive solution 
algorithm to find the optimal harmonic coordinates that define 
the LS solution for our Wavelet waveform. 
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FIG- 5 io the Matlab 5.0 code for the LS algorithm used to 
design the Wavelet waveform in FIG. 6. 



5 

FIG. 6 plots the dc PSD in dB units for the frequency 
response of ^he — our m other Wavelet at dc and the square-root 
raised-cosine (sq-rt r-c) waveforms with excess bandwidth 
10 parameter a= 0.22, 0.40 vs. the normalized frequency in units 
of symbol rate 1/T S . 

FIG. 7 plots the dc PSD in dB units for our waveforms vs. 
the normalized frequency fT s for the our new W avelet waveform 
15 designed with modifications to the LS algorithms for application 
to constant amplitude BEM communications, and for a Gaussian 
minimum shift keying (GMSK) waveform. 

FIG . — 1 — ploto — the — de — P&B — ana — dB — unito — £e* — the — frequency 
20 response of the mother Wavelet at dc and an optimized Gaussian 

minimum — shift — keying (GMSK) waveform ve-s tfee — normalized 

frequency f T a — 

FIG. 8 plots the amplitude of the dc radar ambiguity 
25 function vs. the normalized frequency fT p and normalized time tT c 
where the pulse time T p and chip time T c are both equal to the 
communications symbol interval T s for this example, for the new 
waveform and an unweighted chirp waveform. 

30 



35 
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DI S CLO S URE OF INVENTION 

DISCLOSURE OF THE INVENTION 



5 

Wavelet multi resolution waveforms incorporate a frequency 
translation/ — arc complex, — arc designed in frequency harmonic 
coordinates , — a single mother Wavelet design is used for all 
scales, — and in the t-f space arc realized as FIR filters and 
10 waveforms . 

New Thcsc Wavelet waveforms in this invention disclosure are 
generalizations of Wavelets in t-f space which enable them to be 
useful for communications and radar applications. This 

15 generalization is accomplished by 1) the addition introduction 
of a frequency translation, by 2) relaxing — changing the 
orthonormality condition in equation (4) to apply to waveforms 
within the same space {q} and over the scales {p} with the 
inclusion of the frequency translation, and 3) by their 

20 characterization and design in the Fourier domain. With frequency 
translation the analytical formulation of these new waveforms as 
a function of the baseband or mother waveform centered at dc( dc 
refers to the origin f=0 of the frequency space) becomes 

25 

New waveform as a function of dc waveform 

if/ pqr («) = T* 1 y/(T» n - qM ) e* 2 **** 7 ( 7 ) 

30 where f c (P/r) is the center frequency of the frequency translated 
dc waveform, at scale "p" and frequency index "r". The purpose 
of the frequency index "r" is to identify the center frequencies 
of the waveforms at the scale w p" in the t-f space. The dc or 
baseband or mother waveform is — fcke — generalization of — the mother 

35 Wavelet for multi - resolution waveforms. 
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These waveforms satisfy the complex orthonormality 
equations 



Orthogonality equations 

Y\i/ = 1 iff p = kandq = m and r = v 

^ Y p,q,r Y k,m,v 

= 0 otherwise (8) 



where is conjugation, and are generalizations of the 

orthonormality equations for the analytical Wavelets in (4) . 

The new waveforms in equation (7) expand the Wavelet 
analytical formulation to include a frequency variable. 
Wavelets are functions of the scale and translation parameters 
"p , q" . The concept of an additional parameter to provide an added 
degree of flexibility in their tabulation was first introduced by 

feke — Coif man — library — o4 Wavelets — £-2-3 — which — th&e — aft — oocillaton 

parameter £hsh? — tabulation — that roughly corresponds fee the 

frequency of oscillation of the Wavclcto. The frequency variable 

together with the Fourier domain design are entirely new means 
for deriving these new waveforms as generalization of the 
traditional Wavelets a**d — their — modification — with — fc£e — added 
parameter in — [2] . 

Now These waveforms are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar in the t-f space. The basis vectors for 
this metric space V consist of a subset of the admissible set 
of scaled and translated waveforms { y/ pq r , Vp,q,r } derived from the 
dc waveform y/ as per equation (7) . An admissible waveform is 
any combination that covers K=t-f space. We are interested in 
the Fourier domain representation of the dc waveform y/ in V , 
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and in particular in a subset of the discrete Fourier transform 
(DFT) harmonic coefficients over the Fourier domain which we 
intend to use as the design coordinates. Starting with the z- 
transform and continuous Fourier transform, the DFT harmonic 
coefficients are defined by the following equations 



DFT Harmonic Coefficients (9) 
y/(z) = ^i//(n) z' n z-transform 

n 

\f/(co) = ^^(n) e~ l<m Fourier transform 

a 

y/ k = ^y/(n) WjJ" DFT harmonic coefficients for \/k 

n 

where 

k n 

= YjV* sin((6?/2- ^k/N , )AT')/ Nsin(oV2 -^k/N 1 ) 

k 

= ^ y/ k [Harmonic interpolation for M k"] 



k 

i2/ricn/lsr 



{k} = DFT frequency or harmonic coefficients such that f k = k/N'T 

where f k is the harmonic frequency corresponding to "k" 
N f = length of y/(n) 

where ^(co) = \i//(o))f is the power spectral density of the dc 
waveform. These equations define the frequency representation of 
our new waveforms in terms of the available set of harmonic 
coefficients, which set is considerable larger than required for 
most applications . 

New — Our new w aveforms are an orthonormal basis with no 
excess bandwidth, — enee — propcrtioo — that — we — will — dcmonotratc . 
Orthomalitv and no cxccoo bandwidth aro which properties which 
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are asymptoticly approached by our new waveforms to within design 
accuracies inherent in communications and radar. To procood 
demonstrate these properties we need to identify the structure of 
the dc waveform in V . We start with definitions for the 
5 parameters and coordinates in the following equation (10) . The 
waveforms derived from the dc waveform will be designed to be 
orthogonal over both time translates "MT" and frequency 
translates "1/£t¥MT" which respectively correspond to the Wavelet 
symbol spacing T S =MT and the adjacent channel spacing 1/T S . This 

10 means the orthogonal spacing of the waveforms in V are at the 
time-frequency increments (MT, 1/iiTMT) - (T s . 1/T S ) . In the 
interests of constructing our f*ew — orthonormal multi-resolution 
waveforms to cover V it will be convenient to assume that M,L 
are powers of 2. -We need the following definitions for the 

15 parameters and coordinates. 



20 



25 



30 
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Parameters and Coordinates 



N' - Length of \|/ which is an even function about the center and 
which spans an odd number of points or samples 
= ML + 1 where M, L are assumed to be even functions for 

convenience of this analysis 
= Number of points of y 
M = Sampling interval for y/ 

= Spacing of y/ for orthogonality 
L = Length of \|/ in units of the sample interval M 

= Stretching of vf over L sample intervals 
n = n 0 + n,M 

= partitioning into an index n 0 over the sample 

length n 0 = 0,1,.. .,M- land an index n, =0,1...,L-1 
over the sample intervals 
k =k 0 + k,L 

= partitioning into an index k 0 over the harmonic 
frequencies k 0 = 0,1,...,L-1 corresponding to the 
stretching and an index k, = 0,1,...,M over the harmonics 
frequencies corresponding to the admissible frequency 
slots for y/ 

The harmonic design coordinates are selected using the 
following observation. For most applications and in the following 
development it is assumed that the waveforms are spectrally 

10 contained in the frequency interval MT corresponding to the 

frequency spacing. This suggests the harmonic design coordinates 
be restricted to the subset of L harmonics { k 0 =0, 1 , . . . , L-l } 
covering this spacing. These L harmonics correspond to the 
stretching of the mother waveform over the L repetition 

15 intervals. 
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Obviously, for some applications as will be demonstrated 
later, the spectral containment is spread out over several 1/LT 
frequency increments whereupon one must increase the subset of 
5 design harmonics to possibly 2L, 3L or larger. 

Pursuing — this — application — using — L — design — harmonics , The 
•tkeDFT equations for the dc waveform in (9) when rewritten in 
terms of the L harmonic design coordinates Vk 0 } become: 

10 

DFT equations for dc waveform 

y/ k = ^y(n) WjJ° n harmonic design coordinates 

n 

15 ^(n) =(1/N') W^° n new waveform defined in terms (11) 

ko ° 

of the L harmonic design coordinates {^ o ,Vk 0 } 

It will now be shown that It will bo proven that the use of 
20 these L harmonic design coordinates is sufficient -fet — ewe — aew 

waveform under time translates to be a basis for the 

corresponding subspace of V— This — which m eans these waveforms 
provide a complete set of coordinates to describe this subspace. 
We will use the 

25 will use the theorems of Karhunen-Loeve and Mercer and will 
limit the demonstration to the dc waveform for simplicity and 
without loss of generality. We start by considering Consider the 
expansion of a random complex sequence {z(n), Vn} in a series of 
waveform coordinates consisting of time translates of the mother 

30 waveform iff . The sequence {z(n), Vn} is a zero-mean stationary 
random process which is orthonormal over the sample interval "M" 
and has a frequency spectrum which is flat and extends over the 

19 



frequency range 1/LT which is centered at baseband corresponding 
to a zero frequency. This means the {z(n), Vn) cover the 
subspace of V corresponding to the scale of our new dc waveform 
y/ and its time translates { y/(n -qM) = y/ q (ri), Vq } . In addition, 

5 the V is now considered to be extended over a time interval 
which is relatively large compared to the N-dimensional t-f space 
in FIG. —1,2 to avoid end-effects on the analysis. We start by 
approximating the sequence {z(n), Vn} by the {z(w),Vn} where: 

10 

z(n) = SZ q Kn-qM) 

- ZZ q ^n 0 + (^-tfM) 

where the complex coefficien ts (z q }are 

derived from the original sequence using the ( 12 ) 

orthogonality properties of y/ 
Z q =Iz(n) ^*(n-qM) 

n 

= £z(n)^ q (n) 

D 

The following equations prove that the coefficients {Z q , Vq} are 
orthonormal : 

15 

Z q Z q . * = Iz(An-qM)z*(An-q , M)^ q (An)^ q . *(An) (13) 

An 

= ^ qq ' E|^ q ( n )| 2 since the sequence {z} is orthonormal 

n 

for qM time translates 
= 5. with normalization of the energy of y/ 



20 Equations (12) and (13) together prove the Karhunen- 

Loeve's theorem which proves the following equation for the 
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accuracy in approximating the stochastic sequence {z(n), Vn} by 
{z(w),Vn}. This accuracy is expressed by the expected "E(o)" 
squared error w (o)" in this approximation: 



We need to prove that the right hand side of this equation is 
zero which then proves that the approximating sequence is equal 
10 to the original sequence in the mean-square sense. In turn this 



the original sequence {z(n), Vn} which is our goal. 

The right hand side of equation (14) when set equal to zero 
15 expresses Mercer's theorem so our goal is to prove Mercer's 
theorem. To do this we use the DFT of if/ in equation (11) and the 
coordinates in (10) to evaluate the right hand side of equation 
(14) . We find 



5 




(14) 



proves that the new waveform coordinates { if/ g , Vq } are a basis for 



i - + «,M)| 2 = l-Jj;^, w^r 




where T = (VL^Wf^ 
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= sin(;rAk 0 )/Lsin(;r Ak 0 /L) 
= 1 for Ak 0 = 0 



(15) 



= 0 otherwise 
This proves that 



1- +n t M)\ =0 Vn 0 



which proves that equation (4*14) reduces to 
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E[z(n)-z(n)\ 2 \= 0 



5 which as per the above proves that the set of multi-resolution 
waveforms is a basis. This proof easily generalizes to the 
multi-resolution waveforms at all of the scales {p} and time 
translates {q}, and to the expansion of the harmonic design 
coordinates over 2L, 3L, ... as required by the application. 

10 

New Next we demonstrate that these waveforms have multi- 
scale properties^ that we will demonstrate. We start with the 
observation that the ¥i*e — Fourier domain design for the new 
waveforms provides a natural and easy way to derive the complete 
15 set of waveforms { y/ pq , Vp,q } for the space V from the design of 
the dc waveform y/ by using the same invariant set of Fourier 
domain harmonic design coordinates {NE^ ,Vk 0 } derived for the dc 

multi-resolution waveform. This demonstration requires that we 
show 1) how the multi-scale transformations are implemented with 
20 the design in the Fourier domain, and 2) how the waveform design 
remains invariant under scale changes. 

First consider the multi-scale transformation which derives 
the waveforms at the scale and shift parameters "p,q" for from 

25 the dc waveform at scale "p=0" and centered at the origin "q=0". 
We begin by extending the parameters and coordinates in equation 
(10) , to include both scaling and subsampling or decimation in a 
form that is equivalent to the iterated filter bank construction 
which is used to derive the current Wavelet w aveform using the 

30 filter scaling functions. Starting with the coordinates at scale 
"p=0" the parameters and coordinates at scale "p=p" are given by 
the equations: 
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^£=0 parameters and coordinates 

n = n 0 + ni M (17) 
n 0 = a 0 + ai 2 + ... + a m _i 2 m_1 
= M points 
10 M = 2 m 

ni = b 0 + bi 2 + ... +bi-i 2 1 " 1 

= L points spaced at M sample intervals 
L ~ 2 



15 p=p parameters and coordinates 

2 _p n(l2 p ) = scaled by "2 _p " 

and subsampled or decimated by 2 P :1 

= n 0 (p) + ni(p)M 

2" p n Q (l2 p ) = scaled by "2" p " 
20 and subsampled or decimated by 2 P :1 

= n 0 (p) 

= a p + a p .!2 + ... + ap^-i 2 p+m - x 

= M points spaced at 2 P sample intervals 

2 _p ni(i2 p ) = scaled by "2' p " 
25 and subsampled or decimated by 2 P :1 

= ni(p) 

= bp + b p+1 2 + ... + b p+1 -! 2 P+W 

= L points spaced at M2 P sample intervals 



30 

together with the observation that the sampling interval "T" is 
increased to "2 P T" under the scale change from "p=0" to "p=p" and 
subsampling or decimation from "1:1" to "2 P :1". Combining these 
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equations with the analytical formulation in (7) and the Fourier 
domain representation in (11) enables the waveforms at the 
parameters "p,q" to be written as a function of the Fourier 
domain harmonic design coordinates: 

5 

V„M =2- lV2 ^(2- p n-qM)e i2 ^ r)nT 



I r / \ \a\ ilnf e {p,r)n(p)2*T (18) 



10 for all admissible scale, translation, and frequency index 
parameters "p, q, r" . 

Next we need to demonstrate that the frequency domain 
design in (11) remains invariant for all parameter changes and in 

15 particular for all scale changes. This multi-scale property 
expresses the accordion behavior of the design in that the 
Wavelets at different scales are simply the stretched and 
compressed versions of the mother waveform with the appropriate 
frequency translation indices. This multi-scale invariancy means 

20 that the design for a M=16 channel filter bank remains the same 
for M=100 or M=10,000 channel filter banks, when the overlap L 
and the performance goals remain constant. To demonstrate this 
invariant property across scales, we consider the hh? — multi- 
resolution waveform at scale "p" with the other parameters set 

25 equal to zero for convenience "q=0, r=0" and without loss of 
generality. The Fourier domain frequency response can be 

evaluated starting with the original formulation in equation (7) : 



30 
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DFT at 



q=0, r=0' 



(19) 



= 0>N*) 5X Y W n- 

k 0 n(p) 



-{JN'2 p T -k 0 )n(p) 



k 0 



sin(^r (fNT2 p T-k 0 )) 
N'sin(;r (fN'2 p T-k 0 )/N ! ) 



= J]*F^ [Harmonic interpolation for "k 0 "] 



10 

This only differs from the harmonic representation in equation 
(9) in the restriction of the design coordinates to the subset of 

harmonic coefficients {^,Vk 0 } and the stretching of the time 

interval to "2 P T" corresponding to the scale "p=p". The harmonic 

15 interpolation functions are observed to remain invariant over 
scale changes upon observing that the frequency scales as 
"f~l/2 p T" which means the "frequency* time" product remains 
invariant with scale changes as per the fundamental property of 
the waveforms. This means the harmonic interpolation functions 

20 remain invariant with scale change and therefore the frequency 
response remains an invariant. This demonstrates the waveform 
design is an invariant across the waveform scales which means we 
only need a single design for all scales or resolutions of 
interest. 

25 

LS design algorithms for new waveform will be described to 
illustrate the advantages our ftew — waveform has over current 
designs. The two LS algorithms described are the eigenvalue and 
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the gradient search which respectively can be reduced to 

algorithms which are equivalent to ^he original current 

eigenvalue -FH — and Remez-exchange -{-8-1 — waveform design algorithms 
for application to a uniform filter bank. We consider the t-f 
5 space which is spanned by a uniform polyphase filter bank 
consisting of M channels at the frequency spacing f w =l/MT where T 
is the digital sampling interval, and the filter waveform FIR 
time response is stretched over L sampling time intervals T s . 
This polyphase filter bank is ideally decimated which means the 
10 filter ouput sample rate 1/T S is equal to the channel-to-channel 
spacing T S =MT, equivalent to stating that there is no excess 

bandwidth a=0. Our design for this topology is immediately 
applicable to an arbitrary set of multi-resolution filters 
through the scaling equation (18) which gives the design of our 
15 waveform at arbitrary scales in terms of our design of the dc 
waveform. 

For this polyphase filter bank used to construct the dc 
waveform or filter impulse response, our LS example design 

20 algorithms will use 5 metrics consisting of the 2 prior art 
passband and stopand metrics, and the 3 new metrics consisting 
of the ISI, ACI, and QMF, and solve the LS minimization problem 
using as design coordinates the subset of harmonic coordinates 
which are a basis. Since our 2 — two example LS design algorithms 

25 only differ in the use of an eigenvalue LS optimization and the 
use of a gradient search LS optimization, the flow diagrams for 
the construction of the cost functions and the solution for the 
optimal waveform will be identical. However, there are 
differences in the construction of the cost functions from the 

30 respective metrics and in the iterative solution mathematics. 
Both LS solutions for the harmonic design coordinates minimize 
the weighted sum of the error residuals or cost functions from 
the 5 metrics. These design coordinates fe* — are the Fourier 

harmonics {^,Vk 0 } for {k 0 = 0, 1,...,L-1} . Resulting algorithms 
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are easily extended to the applications requiring the design 
coordinates to cover 2L, 3L, ... harmonics. 



Frequency domain design coordinates are related to the 
5 waveform time domain digital samples or coordinates as follows. 

Mappings of time Q — frequency to time (20) 

Time domain design coordinates {^(w), Vnjare real and 
10 symmetric and can be represented by the reduced set {h t (n), 

n=0,l,...,ML/2} 

h t (n) = ^(0) for n=0 

= 2y/(n) for n=l, 2, ML/2 

= time domain design coordinates 

15 

Frequency domain harmonic design coordinates {^ o ,Vk 0 } are 
real and symmetric and can be represented by the reduced 
set {h f (k) , k=0,l,,...,L-l} 

h f (k) = tf/ ko for k=k o =0 

20 == 2y/ ko for k=k 0 =l, 2, L-l 

= frequency domain design coordinates 

Mapping of the frequency coordinates {h f (k), k=0 , 1 , L-l } 
into the time coordinates {h t (n), n=0 , 1 , ML/2 } is defined 
25 by the matrix transformation to within a scale factor 

h t = B hf 

where 

h f = (h f (0), h f (ML/2) ) t transpose of column vector 
h t = (h t (0), h t (ML/2) ) t transpose of column vector 
30 B = (ML/2 +1) x L matrix 

= [ B kn ] matrix of row k and column n elements B kn 
B kn - 1 /ML for n=l 

— 2 cos(27tkn/ ML) otherwise 
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Mappings — — fcke — frequency — coordinated — (-h ^(k) > — k~Q ,1 , ...,L 1) 
into — tfee — time coordinates — (-h ^ ( n ) , — n~0, l,...,ML/2) — a«d — vice 
vcroa arc defined by the matrix tranof ormationo 

h t =— B fef 

h^ — = — B — h^ 



where 

fe^— = — fo ^ ( 0 ) , — „. / h^(ML/2) )* — tranopooc of column vector 

10 k^— ~ (hfe(0), ...,hfe(ML/2) )* — tranopooc of column vector 

fi — (ML/2 i 1) x L matrix 

= — [ — B^— | matrix of row 1c and column n elements 



1 /ML + 1 for n-1 

2 coo (2rckn/ ML + 1) othcrwioe 



15 B~ * =• L x (ftL/2 + 1) matrix 

-i 



+Bj^ matrix of row k and column n clcmcnto -B^- 



for n-1 

2 coo(27ikn/ ML + 1) otherwioc 

and 



20 B~^~B - I L x L identify matrix 

B—B" 4 — =^— I (ML/2 I I) ac (ML/2 i l) identity matrix 

wherein the N'=ML+1 has been replaced by ML since a single end 
point has been added to the the FIR to make it symmetrical for 
25 ease of implementation for the example Wavelet being considered 
with a sample at the mid-point that makes the number of samples 
N' an odd number. 



30 Passband and stopband metrics and cost functions are 

derived with the aid of FIG. 3 which defines the po wer opcctral 

denoity (-PSDi- parameters of interest for the passband and 

stopband of the PSD ^(co) for communications applications. 
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Requirements for radar applications include these listed for 
communications. Referring to FIG. 3 the passband 11 of the 
waveform PSD is centered at dc (f=0) since we are designing the 
dc or baseband waveform, and extends over the frequency range <0p 
5 extending from -co p /2 to +<0p/2 12 in units of the radian 

frequency variable co=27tfT 13 where T io the digital sampling 

interval defined in FIG. — The frequency space extends over the 
range of f=-l/2T to f=+l/2T which is the frequency range in 
FIG 1 translated by -1/2T so that the dc waveform is at the 

10 center of the frequency band. Quality of the PSD over the 
passband is expressed by the passband ripple 14. Stopband 15 
starts at the edge 16 of the passbands of the adjacent channels 
+ /-co a /2 16 and extends to the edge of the frequency band co=+/- 
7i 17 respectively. Stopband attenuation 18 at +/-co a /2 

15 measures the PSD isolation between the edge of the passband for 
the dc waveform and the start of the passband for the adjacent 
channels centered at +/-co s 19. Rolloff 20 of the stopband is 
required to mitigate the spillover of the channels other than the 
adjacent channels, onto the dc channel. Deadband or transition 

20 band 21 is the interval between the passbands of contiguous 
channels, and is illustrated in FIG. 3 by the interval from G)p/2 
to ca a /2 between the dc channel and adjacent channel at (o a . 
Waveform sample rate (0 3 22 is the waveform repetition rate. 
For the LS example algorithms, the waveform sample rate is equal 

25 to the channel-to-channel spacing for zero excess bandwidth. 
Therefore, 1/T S = © s /27tT=l/MT which can be solved to give Q) s =27i/M 
for the radian frequency sampling rate of the filter bank which 
is identical to the waveform repetition rate. 

30 We start by rewriting the DFT equations for the dc waveform 

in (11) as a function of the {h f (k), k=0, 1, L-l } 
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y/{co) = ^^(n)cos(no>) 



n 



= c* B h f using (2 1 ) and the definition of the vector V 



(21) 



c = ( 1 , cos(fi)), . cos( (ML/2) 6)) y transpose of column vector 

which is equivalent to the equation for ys(p)) in (9) expressed as 
5 a linear function of the {h f (k), k=0, 1,...,L-1} . However, — thio 
functional — form ouggcotcd by the — eigenvalue — formulation — R-) — is- 

more — convenient — te — analyze . An ideal tt c" vector "c r " will be 

introduced for the passband and the stopband in FIG. 3— in order 

to identify the error residual Sy/((o) at the frequency "o?" in 
10 meeting the ideal passband and stopband requirements. The ideal 
PSD is flat and equal to "1" for the passband, and equal to "0" 
for the stopband. We find 

15 Error residuals for passband and stopband (22) 



The LS metric for the passband and stopband can now be 
constructed as follows for the eigenvalue and the LS optimization 
25 (or equivalently, the LS algorithm) design algorithms— 



c r = (l f l,... r l) t passband ideal xx c' 
(0,0 f ...,0) t stopband ideal "c 
Sc = c r - c error vector 
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Sy/((D) = & l Bh f 

= Residual error in meeting the ideal spectrum at 'V 



.it 



30 



30 



Passband and stopband metrics 



(23) 



J (band) = * dfi> Eigenvalue 



band 



" IK i ? 



Eigenvalue 
LS 



(23) — continued 



—where 



IN - 



passband 
stopband 

— |(B t && t B)dfi> 

band 

= L x L matrix 

= ( 5\|/(<0i) 8\|/(a>u) ) t 
= vector of error residuals at the 

frequencies a>i , ©u across the band 

norm or length of the vector (o) and which 

includes a cost function for the errors of 
the individual components 



band = [0, (Dp ) 

= (®s f *] 

1 

R - 

band 

5y 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 
coordinates {h f (k), k=0, 1, L-l } whereas with the LS approach it 
is sufficient to give the LS metrics as vector norms with 
imbedded cost functions or an equivalent formulation. 

QMF metrics express the requirements on the deadband that 
the PSD's from the contiguous channels in FIG. 3 add to unity 
across the deadband [co p , a> s ] in order that the filters be QMF 
filters. By suitable modification of the error vector 8c, the 
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previous construction of the passbanci and stopband metrics can be 
modified to apply to the deadband. This is a Wavelet rcquiromcnt 
eft — the — frequency — coverage — ef — the — Wavelet — baoio — ets — observed — if* 
Table 1. We find 



Deadband metrics (24) 

10 



J (deadband) = (l^(^y)i do) Eigenvalue 

deadband d J aid 

(24) continued 
= h' f R h f Eigenvalue 

= ||<^ I 1 LS 

15 where 

8c = c r - c(©) - c(7i/M - (o) 



where c(co)=c as defined in (22) and (23), and c(7t/M - co) = c at 
20 the offset frequency x 'ti/M - to" corresponding to the overlap of 
the contiguous filters over the deadband . 

Orthonormality metrics measure how close we are able to 
designing the set of waveforms to be orthonormal over the t-f 

25 space, with the closeness given by the ISI and the ACI . ISI io 
the — non- orthogonality — error — between — channel — output — samples 

separated — by — multiples — e£ — the — sampling — interval 1/MT — seconds 

where — T — ±-s — the — sample — time and M is — the — interval — e£ — contiguous 
samples . — AGi — is — the — non - orthogonality — error — between — between 

30 channel — output — samples — within — a — channel — aftd — the — samples — i-fr 
adjacent — channels — at — the — same — sample — time — aftd — at — sample — times 
separated by multiples — e-f — the — sample — interval . As — observed as 
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noise — contributions within each sample — in a given channel, feke 

161 is the noioc contribution due to the other received waveforms 
at the different timing offsets corresponding to multiples of the 

sampling interval. Likewise, the AC1 is the noise contribution 

5 €tee — — the — other — waveforms — i** — adjacent — channels — a^fc — fc£e — same 
sampling time and at multiples of the sampling interval. 
ISI and ACI errors are fundamentally caused by different 
mechanisms and therefore have separate metrics and weights to 
specify their relative importance to the overall sum of the LS 

10 metrics. ISI is a measure of the non-orthogonality between the 
stream of waveforms within a channel as per the construction in 
FIG. 3. On the other hand, ACI is a measure of the non- 
orthogonality between the waveform within a channel and the other 
waveforms in adjacent channels. This means the stopband 

15 performance metric has a significant impact on the ACI due to the 
sharp rolloff in frequency of the adjacent channel, and the ACI 
metric is then a measure of the residual non-orthogonality due to 
the inability of the stopband rolloff in frequency from 
completely eliminating the ACI errors. 

20 

We assume that the received waveform is identical to the 
filter waveforms and is transmitted at the filter output sample 
intervals equal to -t/MT seconds. The second assumption means we 
are assuming the receiver is synchronized with the received 
25 signal. Since there is no information lost by sampling 
asynchronously with the received waveform, we are free to make 
this synchronization assumption without loss of generality. ISI 
metrics are derived in the following set of equations. 

30 ISI metrics (25) 

Mapping of h f into \|/ 

vj; = (v|/ (-ML/2) , ... , \|/(ML/2)) C transpose of column vector 
= H h f 
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H 

Hkn 



= (ML+1) x L matrix of elements H kn 

= 1 for n=l 

= 0.5 cos (2nkn/ML+l) otherwise 



5 Offset matrix A 

A = L x (ML+1) matrix of elements A kn 

A kn = 0 for row vector k=l 

= [0 0 ... \|/(-ML/2) ... y( (L-k)M+l) ] row vector k 
10 (25) continued 

12 kM LM+1 column elements 

ISI error vector 8E 
5E = L x 1 column vector 
15 == A H h f 

ISI metric 

J (ISI) = 5E fc 8E Eigenvalue 

= Non-linear quadratic function of h f 

20 = WseU ls 
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ACI metrics are derived using the ISI metric equations with the 
following modifications. 

ACI metrics (26) 



Mapping of h f into i|/ is the same as developed for ISI 
Offset matrix A elements are changed as follows to apply to 
30 Channel 1: 

A kn = 0 for row vector k=l 

= [0 0 ... \|/(-ML/2)0^ ... v((L-k)M+l) W^ L ' k) ] 
1 2 kM+1 LM+1 



34 



which means the ACI error vector 8E is 
5E = L x 1 column vector 

= AH h 4 h f 

5 ACI metric for the two contiguous channels 

J(ISI) = 2 8E t 8E Eigenvalue 

= Non-linear quadratic function of h f 

= 2 ||<?£||- LS 

10 where the factor "2" takes into account there are two contiguous 
channels or one on either side of the reference channel 0 in FIG. 
3. Because of the fast rolloff of the frequency spectrum the 
addition of more channels into the ACI metric is not considered 
necessary, although the functional form of the ACI metric in (26) 

15 allows an obvious extension to any number of adjacent channels 
which could contribute to the ACI. 

Cost function J for the LS algorithms is the weighted sum 
of the LS metrics derived in (23) , (24) , (25) , (26) . The LS 

20 algorithms minimize J by selecting the optimal set of frequency 
coordinates {h f (Jc),Vk} for the selected set of parameters used 
to specify the characteristics of the dc waveform, frequency 
design coordinates, LS metrics, and weights. Cost function and 
optimization techniques are given by the equations 

25 

Cost function J (27) 
j = ^w(metrics) J(metrics) 

X 

= weighted sum of the LS metrics J (metrics) 

30 where 

metrics = passband, stopband, deadband, ISI, ACI 
{w (metrics) , ^/metrics} = set of weights 
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^ w(metrics) 

JC 



= 1 normalization 



Optimization goal 

Goal: minimize J with respect to the selection of 



5 



the {h f (k) , Vk} 



Optimization algorithms 

Two algorithms are the Eigenvalue and the LS optimization 
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where the eigenvalue optimization algorithm uses the non- 
linear quadratic formulations of the LS metrics and the LS 
optimization algorithm uses the norm formulations for the 



LS metrics. 
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FIG. 4 is a summary of the LS metrics and the construction 
of the cost function J. Design parameters 23 are the input and 

output — design — parameters . Input — parameters — a^e — tke — number — 

polyphase — channels — M — e* — cquivalcntly — fefee — number — e# — digital 

20 samples at spacing T over the symbol interval T g ~MT , the length 

of the — FIR time — response — for the waveform in units — e£— L,_ which 
a^e — the — number — e-f — digital — samples — pe^e — wavcf orm — repetition 
interval T s so that the total number of digital — samples — for the 
symmetric FIR time response is equal to N'= ML+1, number of DFT 

25 samples per FIR length n_fft for implementation of the LS 
algorithms, passband radian frequency (Dp, stopband radiam 
frequency (i) ar waveform repetition rate in radian frequency co s/ 
selection of the set of design coordinates {h f } to be used in 
the optimization, and the metric weights {w (metrics) } . Output 

30 parameters are the set of harmonic design coordinates {h f } that 



minimize J. Band metrics 24 are the passband, stopband, and 
deadband metrics defined in equations (23) , (23) , (24) 
respectively. Interference metrics 25 are the ISI and ACI 
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metrics defined in equations (25) and (26) respectively. LS cost 
function J 26 is the weighted linear sum of the metrics defined 
for the band 24 and the interference 25 as defined in 
equation (27) . 

5 

FIG. 5 is a flow diagram of the LS recursive solution 
algorithm. There are two loops with the topology constructed so 
that the outer loop 27 is an iteration over the set of metric 
weights {w (metrics) } , and the inner nested loop 28 is the 
10 recursive or iteratived LS solution to find the optimal {h f } for 
the given design parameters and weights where each step refers to 
an iteration step in the solution. A recursive LS solution is 
required to be able to solve the highly non-linear interference 
metrics as well as the band metrics when one chooses to use non- 
15 linear techniques to construct these metrics. 

Consider the inner nested loop 28. The recursive LS 
solution starts with the initial step i=0 29 which begins with 
the selection of the input design parameters and weights and the 

20 selection of an initial set of values for the {h f } 30. Next the 
band metrics are calculated 31. Since this is the initial step 
i=0 32 the cost function J is restricted to the linear band 
metrics in order to fine the linear approximation to {h f } 33 to 
initialize the non-linear solution starting with step i=l 34 

25 wherein the highly non-linear interference metrics are included 
in J. Following the signal flow, for step i=l the band metrics 
31 and interference metrics 35 are calculated, weighted, and 
summed to form the cost function J 36. An iterative solution 
algorithm 37 finds the best approximation in step i=l to the 

30 {h f } which minimizes J using the approximation to {h f } from the 
previous step i=0 to linearize the search algorithm coefficients 
for the eigenvalue and LS step i=l iteration. If there is no 
convergence to the correct solution for {h f } this recursive 
sequence of calculations is repeated for step i-2 34. This 

35 recursive solution technique is repeated for subsequent steps 

37 



until there is convergence 38 whereupon one exits this inner 



loop . 

Consider the outer loop. After exiting the inner loop 38 
5 the solution for {h f } is tested to see if it meets the 
performance goals 39. If not, a new set of metric weights is 
selected and the inner loop is initialized i=0 40 and the 
inner loop is used to find the next solution for {h f } . This 
process continues until the performance goals are met or 
10 adequated approximated 39 whereupon the algorithm is exited 
with the final solution set {h f } 41. 

g£€^ 5A,5B,. . .,,arc the Matlab 5.0 oodc for the ij& 
re cursive eigenvalue solution algorithm uocd to design the 

15 Wavelet Wavelet in FIG. 6 . Included arc a listing of the froq eeftey 
domain design harmonics, time response, and data plot ^ftg 
se- ftwarc uocd in the optimization, and an example of the scaling 
of this mother Wavelet into a Wavelet with new oca jrjbftg 
-fd ilation) , time translation, frequency tranolation paramct e^s-r 

20 and the functions used in the code. The Matlab code uoco- aft 
it erative approach to find the optimal frequency domain 
ei genvector solution to the quadratic form for the LS error whose 
e^ror matrix io the weighed sum of the error matrices f or-^fc-ke 
^te pband, paosband, dcadband metrics and the error matrices f ar 

25 the ISI and ACI metrics which arc functions of the eigenvectors . - 
ge r each iteration the ISI and ACI error matrices arc conotruct -ed 
wjr th the frequency domain eigenvectors from the previous - 
j uration. The iteration finds the new frequency doma in* 
e igenvectors which arc then uocd to update the ISI and ACI 

30 mctrico for the next iteration. Thio iteration io repeated unfe jbjfc- 
there is convergence. The Matlab code allows the opcrato ^-^e 
update the metric weights as well as the scenario parameter s— ^ 
#jrnd the optimal choice for the performance goa jr^r 
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FIG 5A Matlab code starts with the selection of the dc s-i-gft 



parameters in 29 Stop 1. Stop 1,1 lists the scenario paramc -fce^ 
wkjr ch arc M~16 digital samples per Wavelet \|/ pomplG intcrval -r 
L~16 nominal length of Wavelet y in units of M, N-ML * 1 Wavelet \\f 
5 length which io represented by N' in equations (10), fo~l io t he 
normalized channel spacing equal to the Wavelet i|/ sample r- a^er 
fp-Q. 88 64 io the paooband frequency interval, n f~16 io the 
number of dcoign frequency harmonics, n fft~102 4 io the number of 
digital samples uood for calculation of the frequency spectrum 

10 centered at 0 frequency. cbno=fib/tte~6 . 0 dB is the design vala e— e£ 
£hc energy per bit E^_to noise power density No-fa ^o, and 
*e- imbal aci-6.0 dB is the assumed channel - to channel pow e^ 
•im balance used in the calculation of the ACI errors. Step 1.2 
de rives the software parameters from the design parameters in 

15 Step 1.1. Step 1.3 lists the optimization paramctcro which a *e 
£k e number of iterations n iteration - 10, and the metric wei ghs 
w- frictric) from equation (27) and in FIG. 4 equal to w(po 9^H= 
e^ pha 1-1. c-2 for paosband, w (stop) -alpha 2-Q.80 for otopb affrdr- 
w(ISI) -alpha 3~2.c 3 for ISI, w (ACI) -alpha 4-Q.5 for ACI, and 

20 w(dcad) ^alpha 5~0 for dcadband. The number of iterations was 
selected to provide a convergent solution for the weighted LS 
ef for metrics w (metric) J (metric) in FIG. 5K 3 8 Step 10 figure (g f 
and for their sum J in figure (1) . Weight values were selected— to 
ep timizc the Wavelet \y filter performance in FIG. 6 and in 

25 figure (3) in FIG. 5L 39 Step 11, Wavelet y ripple, IEI, -AGi 
sjr gnal - to noise SNR power ratio looses in figure ( 4 ), and Wavelet 
y- time rcoponso in figure (S f-=- 

•F^G-= 5B 30 Step 2 arc the initialization calculations prior - 
30 to the start of the iteration loop. Stop 2.1 calculates the on e— 
eir ded Wavelet length parameter m-12 8 . Step 2.2 calculates the 
ma- trix transformation bw matrix which maps the y frequency domain 
de sign eigenvectors into the one - sided time response. Step -2^3 
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rcfcro to the function pmn in FIG. 5P,5Q 4 1 in Step 13.1 which 



ee mputco the LS error matrix which io the sampled data integral 
£j§ c » 8c' ]da> over the band in equation (23) for the the paoobaR d 
ettd otopband metrics J (paoo) , J (stop) . The function pmn -d 
5 computco the corresponding integral over the dcadband in equatio n 
-f24) and io not listed since it wao not used for FIG. 6 . -S£ep 
2rrA constructs a structural matrix c matrix used in the ISI and 
Ad calculationo . Step 2.5 constructs the sample rate templates 
#e r the paoaband upper edge, sample rate marker, and paoob aftd 
10 lower edge for frequency performance evaluation. 

-F£-€-=- 5C 31 Step 3, uses the function pmn to compute the— IrS 
e^ ror matrices passband and otopband in Step 3.1, 3.2. Dcadb aftd 
LS- error matrix dcadband io oot equal to a null matrix in Step 
15 3.3. Step 3. 4 computco the matrix p total equal to the wei ghted 
sum of these matrices. Step 3.5 maps this LS error ma ^i-x 
p-^ total in the time domain into the LS error matrix pw t in the 
-f requency domain using the matrix tranof ormation bw matrix. - 

20 

FIG. 50,50 32 Step 4 initializes matrices used in — the 
j uration loop and otarto the iteration loop for the iterations 
i ^.tcration-1, 10 . Step 4 .1 calculates the eigenvalue — a**d 
ei genvector which minimise the cost function J in equations (27) - 

25 and FIG. 4 whoso LS error matrix is pw t. Step 4.2 uses the 
bw- matrix to map the eigenvector in frequency into the Wav e±e*r 
£4 mc rcoponoc iy-hn and otorco the eigenvector ao the oct — e£ 
Wa velet \y design harmonics ly^-hw cig(k) for 1c-Q, 1, 15 where 
heo been defined in equations (0) and in equations (11) -upeft 

30 recognizing that there arc n f~16~L design harmoni e&^e ^Q^/ - - 
— 15. There is no reason to use the negative harmonics since ^the 
ep-cctrum is real and symmetric about the 0 frequency. Step 4.3 
ee mputco the passband peak to - peak ripple ao xripplc and t **e 
e ^opband attenuation ao xotop ^- 



40 



FIG. 5E 33 Step 5 calculatco the weighted pasoband orre ^ 



metric w (paoo) J (pass) -beta 


paoo, — weighted — stopband — error metric 




and — weighted — dcadband — error — metric 


w (stop) J (stop) -beta stop, — 










format — ie^e — hft — a-ftd — the — LS- 


error matrices pasoband, — otopband, 







dcadband. 



F£€^- 5E,5F 3 4 Stop 6 calculatco the ISI and ACI LS crr ef 
10 matrices, their metrics, and their error contributationo to the 
SN R loos. Step 6.1 calculatco the ISI LS error matrix w matr 4^ 
and the corresponding ISI metric J(ISI)-crrM ioi with Matlab -eede 
wh ich implements cquationo (25) , and calculatco the ISI rcoi -duajr 
eir ror errV ioi that cauoco an ISI SNR loos. Step 6 .2 calcul ates- 
15 the ACI LS error matrix w f matrix and the corrcoponding ACI - 
metric J (ACI)- crrM aci with Matlab code which implcm eft-fe-s- 
e^ uationo (2 6 ) , and calculatco the ACI rcoidual error crrV — aei- 
t& at cauoco a ACI SNR looo. 

20 FIG, 5G 35 Step 7 calculatco the weighted LS error me tres- 



£er — Ij&i — and ACI 


and — the — LS — error matrix — pw t — #e* — the — next 




.1 calculatco the weighted ISI LS error metric 


iteration. — Stop 7 






isi — a-ftd — the — weighed — AGI — LS — error — metric 


w(ISI) J(ISI) -beta 






aci — irft — equationo {2^) — a**d — the — &tm — e4 — the 


w (ACI) J (ACI) -beta 






metrics J~crrM isi. Step 7.2 saves the weighted 


weighted LS error 




ij£ — error — metrics 


and their sum J for each iteration. Step 7.3 







updates the pw t matrix for the next iteration. 

HC. 5C,5H 3 6 Step 8 calculatco the SNR losocs, cn te^es- 



theoc 


looocs — £er — each — iteration, — a-ftd — complctco — the 


iteration 






— pasoband 




— Step — 8-rl — computeo — the — SNR — looo — if* — dB — due — te 








to ACI ao 


ri PP lG 


ao xlooo ripple, — due to ISI as xloss ioi, — due 








8.2 oaves 


xlooo 


aci, — and the total SNR looo ao xlooo total. — Stop 






SNR looses — for each iteration and completes the 


iteration 


thcoo 







35 loop. 
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FSG^ — 5H,5I,5J,5K 37 Step 9 documents the Wavelet 
fee quoncy design harmonics and time response generated by -^fche 
Matlab code for the Wavelet in FIG . \\t 6 . Step Q.l lioto the 
5 Wavelet \y frequency domain harmonic dcoign coordinates 
j fa-tfe ) -hw oig(k) for harmonico k-Q, 1, 15 . Only the positi ve 
harmonics arc listed since the frequency spectrum is rcal ^a-ftd 
sym metric about k-Q . Step 9,2 lioto the Wavelet i|/ time respons e 
\j/(n)-hn(n) for n-Q,l,2,. . , ,128 with the coordinate f-eame 
10 centered at the Wavelet peak. Only the pooitivc digital sampl er 
a^e listed oinco the time response is real and symmetric abou t: 
n=0. 

^€^- 5K 38 Step 10 plots the total weighted LS error mctri e 
15 J in figure (1) and the individual weighted LS error mctrie s— jrn 
£jr ^urc(2) as functions of the iteration number. After about — fr4 
-alt erations both the total and the individual weighted LS crr af 
metrics arc observed to stabilize indicating that the sele cted 
fHHubcr of iterations n itcration-10 io sufficient — 
20 convergence . 

FI&t— 5L 39 Step 11 plots the Wavelet \\f frequency respons e 
irft- figure (3) which is the same as the frequency response in FI G-r- 
G-r- plots the SNR looses due to paooband ripple, ISI, ACI and t he 
25 total loss in figure (4), and plots the time response in figure 

£e r use in the optimization of the LS metric weights in Step 1 ^-3 
ae well as the scenario parameters in Step 3 r-r4r^- 

FIG. 5M 40 Step 12 documents the Matlab code to re scale—the 
30 mother Wavelet \\f in FIG. 6 which was derived in the prcv jrOteg- 
Matlab code in Steps 1 - 11, for now scale (dilation) p, time 
tr anslation q, and frequency translation k parameters compar ed-^to 
^th c values p-Q,q-Q,k-0 for the mother Wavelet \y. An example — &efe 
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o£- new paramctoro documcntcci in the cocio io p~2, q~2 , k~3 and 
with the aooumption that the digital sampling rate remains f-jr^ed 
ffldcr the parameter change. A conotant digital sample rate keeps 
£k e frequency band constant under the parameter chaft§er 

5 

FIG. 4,5 5M 40 Step 12.1 derive the new multi-resolution 
Wavelet \y at baseband. Scaling equations defines the new Wavelet 

^^^ ^ in equation s — (7) , (18) can be used to derive the 

corresponding digital Wavelet at scales p,q,k with parameters 
10 M,L wherein the digital frequency parameter "k" replaces the 
frequency index "r". Equation (7) can be rewritten to define the 
scaling equation which derives this Wavelet \|/ p , q , M .L.k from y: 
res ulting from these parameter changco, in equations -j2Q± 

15 New Wavelet scaled from mother Wavelet (28) 

vR P ,q / M,L ( k_^ 7 ^-z_^,q,k with parameters M,L 

= 2" p/2 y(2" p n-qM) exp(j27i k 2" p n /ML) 

20 There are at least two methods of scaling used to implement 

equation (28) . In method 1 corresponding to the standard Wavelet 
assumption, the subsampling or decimation of n by the factor 2 A p 
stretches the time between contiguous samples by the factor 2 A p 
while keeping constant the number of samples M per Wavelet y 

25 sample interval, and therefore reduces the frequency band by the 
factor 2 A p. For a uniform filter bank, this means the filter 
frequency bands are reduced by the factor 2 A p. In method 2, the 
sampling is held constant corresponding to a constant frequency 
band assumption. The factor 2 A (-p) is divided out in the Wavelet 

30 expression which is equivalent to the M being scaled to 2 A pM 
corresponding to increasing the number of samples per Wavelet 
interval by the factor 2 A p. For a uniform filter bank, this 
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number of channels by the factor 2 A p over the same frequency 



band . 



5 

Case 1 Scaled (dilated) sampling -^ft-w^h fixcd -^4 

oubsamplcs (decimates) n by the factor — 2? 

— ■ Rcproocnt n as 

n - n 0 + 2 P — n— p 
10 n 0 - 0,1, . . 2*-± 

n p , - 0 , 1 /1/ 1 /2, . . . 

■ We observe — S-^fi— p 

* This enables the new Wavelet to be written 

Vp hh* -^ 2--^- n/(n p qM) cxp(j27t k n p /Mi^ 

15 

Ga^ e 2 Sampling n remains fixed and M ^-9- 

— repealed by the facto jE^-^ 

— * New M is M new 
M new - 2^M 

20 ■ This enables the new Wavelet to be wrj r^beft 

\\f^ r ^j==--2 r ^ \|/(n - qM new) cxp(j2ft kn/M ncw L-f 



25 wherein the new Wavelet equation from (7) , (1 8 ) is rc-writfe eft 
using the frequency translation parameter k-r and the cxp jrireir^- 
^ equency translation expression cxp(j2n k 2~ * *n /ML) whie h-^rs- 
re cognized as the discrete Fourier transform kcrna -^ 

30 

In Case 1 corresponding to the standard Wavelet assump ^-j-oftr- 
£k e suboampling or decimation of n by the factor 2 A p strc ^tehes- 
the time between contiguous samples by the factor 2 A p whi 4^ 
4ee eping constant the number of samples M per Wavelet \\f sample 
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jrffrfc crvQl, and therefore rcducco the frequency band by the facto r 
£^-p . The new n can be rcprcocntcd ao the oum of the two fi e4^& 
b— Q — 0,1,. , , (2 A p 1) and n p - 0, 1 / 1, 1 / 2, . . ■ Sc* 3r*ftg- 
e^- cquivalcntly ouboampling rcoulto in the new n equal to n-ftew 
5 - 2 A -p n - n p corrcoponding to a otrotching of the oamplo 
jrft tcrval by 2 A p. The new Wavelet length N new rcmaino cons^ffl^ 
N— rcw-N uoing the n p oampling. 

±n — G aoe 2 which io the communicationo application b eirttg- 
10 conoidcrcd, the oampling io held conotant corrcoponding ^e— a 
ee notant frequency band aooumption. The factor 2 A p io div ided 
ettt in the Wavelet cxprcooion which io equivalent to the M bein gs 
ee alcd to M ncw~2 A p M corrcoponding to incroaoing the numbene— 
oamplco per Wavelet interval by the factor 2 A p. Thio mcano -Gese 
15 2 corrcopondo to incrcaoing the maximum number of channclo by t he 
£e etor 2 A p over the oamc frequency band. The new Wavelet lengt h: 
j^o- equal to N ncw^M ncwL t l -1025 oincc M ncw~2 A pM~64 

•F^G^— 5M 40 Step 12.2 dcrivco the new matrix tranof ormati eft 
20 bw matrix new which mapo the \\f frequency domain de ^-gft 
cigcnvcctoro into the one - oided \y time rcope ft&e-r- 

— 5M,5N 40 Step 12.3 uoco the matrix tranof ormation 
bw- matrix new to map the dooign eigenvector hw cig into the — ftew 

25 Wavelet in three otcpo. In the firot otcp, the mother Wavel et— ±s 
feg- anolatcd into the new baocband Wavcl -et— *t* P , q - o,k - o ^-4* n 0 whie h— 3r& 
£h e mother Wavelet ocalcd by p and with no time and frcqucn ey 
t^ anolationo q-k-Q. The occond otcp tranolatc thio baocba ftd 
Wa velet by qM new to generate -^^ ^q s- hn 1. The laot ot ep 

30 tranolatco thio Wavelet -^^^^^ hn 1 in frequency to give the n ew 
Wa velet ii/ e re7 ^ -=^- hn new. 

Fi€-s— 5N 40 Step 12. 4 ploto the Wavelet time rcoponoco for 
th e mother Wavelet iy-hn and the new Wavelet \\f &r<t ^-~ hn ne w — i-ft 
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£jrg urc(6). Step 12.5 ploto the Wavelet frequency responses fene 
the mother Wavelet n/~hn and the new Wavelet 1^^^^ hn ncw -4ft 
fi gure (7) versus the normalized frequency/ (\y-hn sample rate) — eft^ 
±a— figure (8) versus the normalized f rcqucncy THW r^rr* — ^ hn ne w 
5 sample rate) 

F-^G^- 41 Step 13 documents the functions used in the Matlab 
eod c. Step 13.1 is the code for the function pmn which compute r 
feke- matrix for the J (BAND) in equations (23). FIG. 4 1 Step 
10 lists the code for the function frcq rsp which computes — fcke 
Fo urier transform of the input hn versus the normalized 
^ cquoncy/ (Wavelet \|/~hn sample ratc -H- 

Considcr the inner nested loop 2 8 . The recursive LS 

15 solution starts with the initial stop i~0 29 which begins with 
%h c selection of the input design parameters and weights and the 

se lection of an initial set of values for the {h d 3 0. Ncxt ^fche 

be ftd metrics arc calculated 31. Since this is the initial step 
jr=Q 32 the coot function J is restricted to the linear ban d 

20 metrics in order to fine the linear approximation fe e— R M- 33 to 

in itialize the non-linear solution starting with step i~l 34- 

wh ercin the highly non linear interference metrics arc included 
J. Following the signal flow, for step i~l the band metrics - 
3^- and interference metrics 35 arc calculated, weighted, — a**4 

25 summed to form the coot function J 3 6 . An iterative ooltt feirOft 

algorithm 37 finds the best approximation in step i~l to the 

-fh^ — which minimizes J using the approximation to — f-hd — from the 

previous step i~0 — to linearize the search algorithm coefficients 
for the eigenvalue and LS step i~l iteration. If there is no 

30 convergence — fee — tfee — correct — solution — — R*d — this — recursive 

sequence — e£ — calculations — is — repeated — fe* — stop — 34-= This 

recursive — solution — technique — is — repeated — — subsequent — steps 

until there io convorgonco 38 — whereupon one exits this inner 

loop . 
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Consider the outer loop. After exiting the inner loop — 3S 

the — oolution — fe* (-kg-) — 3rs — tooted — fee — &ee — ir£ — it — mccto — the 

performance goals 39-. If not, — a new got of metric wcighto io 

5 selected — a**d — the — inner — loop — irS — initialized — i=0 4© aftd — the 

inner — loop — irs — uocd — fee — find the — next — oolution — fe* — (-h-rH Thio 

process — continuco — until — the — performance — goals — a-re — met — 

adequated — approximated 39 whereupon the — algorithm io — exited 

with the — final oolution oct — (4*d 4i. 

10 

Applications of this new invention to both communications 
and radar will be given using these example algorithms and other 
algorithms supported by this invention. These new waveforms are 
considered for the applications: l)to replace the square-root 

15 raised cosine waveform (sq-root rc) which is extensively used for 
the third generation (3G) CDMA communications, 2) to replace the 
Gaussian minimum shift keying (GMSK) waveform for constant 
amplitude bandwidth efficient (BEM) applications, and 3) as a 
candidate waveform for synthetic aperture radar (SAR) and real 

20 aperture radar (RAR) applications. 

CDMA communications application for the current and the new 
3G CDMA considers a waveform designed with this new invention as 
a possible replacement for the sq-rt rc waveform with bandwidth 

25 expansion parameter a -0.22 to a = 0.40 . This notation means that 
for a =0.22 the spectral efficiency is (symbol rate/bandwidth) = 
=1/1.22 = 0.82 = 82%. A basic advantage of the waveform 
is the potential for a symbol rate increase within the same 
bandwidth with an increase in the spectral efficiency to « 100% 

30 depending on the application and operational constraints. The dc 
power spectral density or power spectrum (PSD) of the waveform is 
compared to the PSD for the sq-rt r-c in FIG. 6. Plotted are the 
measured PSD in dB 42 versus the frequency offset from dc 
expressed in units of the symbol rate 43 Plotted against the 
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normalized frequency offset are the dc PSD for the new waveform 
44, the sq-rt r-c with a=0.22 45, and the sq-rt r-c with 
<x=0.4 0 46. It is observed that the PSD for the new waveform 
rolls off faster that that for the sq-rt r-c which means that the 
5 new mr waveform will support an increased symbol rate for a given 
available frequency band while satisfying the inherent 
requirements for low ISI and MAI (multiple access interference) . 

Constant amplitude BEM application of the new waveform 
indicates that it is a viable candidate for replacing the current 
preferred modulation waveform which is the GMSK. The GMSK finds 
applications for transmitters which operate their HPA(s) 
amplifiers in a saturation mode in order to maximize their 
radiated power from the HPA(s) , and which require a BEM PSD to 
avoid excessive spreading of the transmitted power. Simulation 
data for the dc PSD is plotted in FIG. 7 for the new waveform BEM 
and the GMSK. Plotted are the measured dc PSD in dB 47 versus 
the frequency offset from dc expressed in units of the bit rate 
48. Plotted against the normalized frequency offset are the dc 
PSD for the new waveform BEM 49 and the GMSK 50, for a length 
parameter L=10 where L is the length of the phase pulse in terms 
of the phase pulse repetition rate. The significance of this 
example data is that the new waveform has the potential to be 
designed to offer a PSD which is less spread out than the current 
GMSK, and therefore an improved BEM waveform. 

Radar RAR and SAR application of the waveform indicates 
that it is a viable candidate to replace the current chirp 
waveforms for wideband signal transmission, when combined with 
30 pseudo-random phase codes. Results of the simulation for the new 
waveform and an unweighted frequency chirp waveform are given in 
FIG. 8. Plotted are the ambiguity function for the new waveform 
51 and the unweighted frequency chirp waveform 52. The dc 2- 
dimensional radar ambiguity function 53 is plotted as a function 
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of the frequency offset in units of fT p and the time offset in 
units of t/T c where T p is the phase-coded radar pulse length or 
length of the phase code and T c is the phase code chip length. 
The chip length is identical to the waveform repetition interval 
5 T s so that T c = T s . It is observed that the new waveform has the 
potential for significant improvements in the ambiguity function 
and by implication in the performance. 

Preferred embodiments in the previous description i-s — are 
10 provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
15 present invention is not intended to be limited to the 
embodiments shown herein fetttand— is not to be accorded the wider 
scope consistent with the principles and novel features disclosed 
herein . 

20 
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FIG. 2 Wavelet Iterated Filter Bank for Tiling t-f 
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FIG. 3 PSD Requirements for Communications 
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FIG. 4 



LS Metrics and Cost Function 
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FIG. 5 LS Recursive Solution Algorithm 

29 



Start 




5 



55 



FIG. 6 PSD for New Waveform and Square-Root Raised-Cosine 




Normalized frequency, symbol rate 43 



5 



10 



15 



56 




58 



APPLICATION NO. 09/826,118 

TITLE OF INVENTION: Wavelet Multi-Resolution Waveforms 
INVENTOR: Orbain A. von der Embse 



Clean version of how the SPECIFICATION 
Will read. 



l 



APPLICATION NO. 09/826,118 

TITLE OF INVENTION: Wavelet Multi-Resolution Waveforms 
INVENTOR: Urbain A. von der Embse 

5 

BACKGROUND OF THE INVENTION 



I. Field of the Invention 

10 

The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data 
communications with data rates up to multiple Tl (1.544 Mbps) and 
higher (>100 Mbps), and to optical CDMA. Applications are 
15 mobile, point-to-point and satellite communication networks. More 
specifically the present invention relates to a new and novel 
approach to the design of waveforms and filters using 
mathematical formulations which generalize the Wavelet concept to 
communications and radar. 

20 



II. Description of the Related Art 

25 

Multi-resolution waveforms used for signaling and/or 
filters and which are addressed in this invention are defined to 
be waveforms of finite extent in time and frequency, with scale 

30 and shift properties of multi-resolution over the time-frequency 
(t-f) space. These waveforms can also be referred to as multi- 
scale waveforms, and include multi-rate filtering and the 
Wavelet as special cases. The emphasis will be on digital 
design and applications with the understanding that these multi- 

35 resolution waveforms are equally applicable to analog design and 
applications . 

1 



Background art consists of the collection of waveform and 
filtering design techniques which can be grouped into six broad 
categories. These categories are: CI) least squares (LS) design 
5 algorithms for filters and waveforms that design to 
specifications on their frequency response, C2) analytic 
filters and waveforms which are specified by a few free design 
parameters that can be sub-categorized into current applications 
and primarily theorectical studies, C3) combinations of CI and 

10 C2 for greater flexibility in meeting communications and radar 
performance goals, C4) special design techniques to control the 
noise levels from intersymbol interference (ISI) and adjacent 
channel interference (ACI) in the presence of timing offsets for 
multiple channel applications, C5) Wavelet filter design using 

15 scaling functions (iterated filter banks) as the set of design 
coordinates or basis functions, C6) filter and waveform design 
techniques for non-linear channels and in particular for 
operation in the non-linear and saturation regions of a high 
power amplifier (HPA) such as a traveling wave tube (TWT) or a 

20 solid state amplifier, and C7) LS dynamic filters derived from 
discrete filtering and tracking algorithms that include adaptive 
equalization for communications, adaptive antenna filters, 
Wiener filters, Kalman filters, and stochastic optimization 
filters . 

25 

Category CI common examples of LS digital filter design are 
the eigenvalue algorithm in "A New Approach to Least-Squares FIR 
Filter Design and Applications Including Nyquist Filters" and the 
Remez-Exchange algorithm in "A Computer Program for Designing 

30 Optimum FIR Linear Phase Filters' . The eigenvalue algorithm is a 
direct LS minimization and the Remez-Exchange can be reformulated 
as an equivalent LS gradient problem through proper choice of the 
cost function- Both LS algorithms use the FIR (finite impulse 
response) digital samples as the set of design coordinates- LS 

35 design metrics are the error residuals in meeting their passband 
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and stopbanci ideal performance as shown in FIG. 3. Category C2 
common examples of analytical waveforms and filters for system 
applications are the analog Chebyshev, Elliptic, Butterworth, and 
the digital raised-cosine, and square-root raised-cosine. For 
5 theorectical studies, common examples are polyphase multirate 
filters, quadrature mirror filters (QMF) , and perfect 

reconstruction filters. Although these theorectical studies have 
yet to yield realizable useful filters for system applications, 
their importance for this invention lies in their identification 

10 and application of ideal performance metrics for filter 
designs. Category C3 common digital example is to start with 

the derivation of a Remez-Exchange FIR filter and then up-sample 
and filter with another bandwidth limiting filter. This results 
in an FIR over the desired frequency band that is larger than 

15 available with the Remez-Exchange algorithm and with sidelobes 
that now drop off with frequency compared to the flat sidelobes 
of the original Remez-Exchange FIR. A category C4 common example 
is to select the free parameters of the category C3 filter in 
order to minimize the signal to noise power ratio of the data 

20 symbol (SNR) losses from the ISI, ACI, and the non-ideal 
demodulation. A second common example is start with a truncated 
pulse whose length is short enough compared to the symbol 
repetition interval, to accommodate the timing offsets without 
significant impact on the ISI and ACI SNR losses. This shortened 

25 pulse can then be shaped in the frequency domain. 

Category C5 Wavelet filter design techniques discussed in 
the next section will serve as a useful reference in the 
disclosure of this invention. Category C6 common example is the 
30 Gaussian minimum shift keying (GMSK) waveform. This is a 
constant amplitude phase encoded Gaussian waveform which has no 
sidelobe re-growth through a non-linear or saturating HPA. 
Category C7 common examples are the adaptive equalization 
filter for communication channels, the adaptive antenna filter, 
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and the Kalman filter for applications including target tracking 
and prediction as well as for equalization and adaptive antennas. 

Minimizing excess bandwith in the waveform and filter 
5 design is a key goal in the application of the C1,...,C6 design 
techniques for communications and radar. Excess bandwith is 
identified as the symbol a in the bandwidth-time product 
BT s =l+a where the two-sided available frequency band is B and the 
symbol repetition interval is T s . Current performance capability 

10 is represented by the use of the square-root raised cosine (sq-rt 
rc) waveform with a= 0.22 to 0.4 as shown in FIG. 6. The goal is 
to design a waveform with a=0 within the performance constraints 
of ISI, ACI, passband, sideband, and passband ripple. This goal 
of eliminating the excess bandwidth corresponds to the symbol 

15 rate equal to the available frequency band 1/T S = B/{l+a) = B for 
a=0. This symbol rate 1/T S =B is well known to be the maximum 
possible rate for which orthogonality between symbols is 
maintained. A fundamental performance characteristic of our new 
waveform designs is the ability to eliminate excess bandwidth for 

20 many applications. 

Scope of this invention will include all of the waveform 
and filter categories with the exception of category C4 special 
design techniques and category C7 LS dynamic filters. Emphasis 

25 will be on the category C5 Wavelets to establish the background 
art since Wavelets are multi-resolution waveforms that can 
eliminate the excess bandwidth and have known design algorithms 
for FIR waveforms and filters. However, they do not have a design 
mechanism that allows direct control of the ISI, ACI, passband, 

30 sideband, and passband ripple. Category C2 theorectical studies 
also eliminate the excess bandwidth. However, they are not multi- 
resolution waveforms and do not have realizable FIR design 
algorithms. So the emphasis in background art will be on 
Wavelets whose relevant properties we briefly review. 
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Wavelet background art relevant to this invention consists 
of the discrete Waveform equations and basic properties, 
application of Wavelets to cover a discrete digital time- 
5 frequency (t-f) signal space, and the design of Wavelets using 
the iterated filter construction. Wavelets are waveforms of 
finite extent in time (t) and frequency (f) over the t-f space, 
with multi-resolution, scaling, and translation properties. 
Wavelets over the analog and digital t-f spaces respectively are 
10 defined by equations (1) and (2) as per Daubechies's "Ten 
Lectures on Wavelets", Philadelphia : SI AM, 1992 



Continuous wavelet 

15 ¥ ( 0 = | fl |-"V(— ) (1) 

a 

Discrete wavelet 

, . | 1-1/2 

a 

where the two index parameters "a,b" are the Wavelet dilation and 
20 translation respectively or equivalently are the scale and shift. 
The \\f is the "mother" wavelet and is a real and symmetric 
localized function in the t-f space used to generate the doubly 

i H /2 

indexed Wavelet \|/ a b- The scale factor "|a| " has been chosen to 

keep the norm of the Wavelet invariant under the parameter change 
25 "a,b". Norm is the square root of the energy of the Wavelet 
response- The Wavelets \j/ a ,b and vj/ are localized functions in the 
t-f space which means that both their time and frequency lengths 
are bounded. The discrete Wavelet has the time w t" replaced by 
the equivalent digital sample number "n" assuming the waveform is 
30 uniformly sampled at W T" second intervals. 
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Wavelets in digital t-f space have an orthogonal basis that 
is obtained by restricting the choice of the parameters "a, fa" to 
the values a=2" p , b=qM2 p where *p,q" are the new scale and 
translation parameters and "M" is the spacing or repetition 
5 interval T S =MT of the Wavelets (which from a communications 
viewpoint are symbols) at the same scale "p" . Wavelets at "p,q" 
are related to the mother Wavelet by the equation 

10 yv q (n) =2" P V(2- p n-qM) (3) 

where the mother Wavelet is a real and even function of the 
sample coordinates. The orthonormality property means that these 
15 Wavelets satisfy the orthogonality equation with a correlation 
value equal to "1". 

XVp,qVk,m =1 iff both p = kandq=m 

n 

= 0 otherwise (4) 

Wavelet representation of a digital t-f space starts with 
20 selecting an N sample time window of a uniform stream of 
digital samples at the rate of 1/T Hz (1/second) equivalent to a 
"T" second sampling interval. The N point or sample t-f space 
in FIG. 1 illustrates a Wavelet representation or "tiling" with 
Wavelets that are designed analytically or by an iterated filter 
25 construction. 

The t-f space in FIG. 1 is partitioned or covered or tiled 

by a set of Wavelet subspaces tyP p9 p = 0,l,...,/w-l} 1 where N=2 m . Each 

Wavelet subspace W p 1 at scale u p" consists of the set of Wavelet 
30 time translations {q = 0,1,..., N/2 P+1 -1} over this subspace. These 
Wavelet subspaces are mutually orthogonal and the Wavelets within 

6 



each subspace are mutually orthogonal with respect to the time 
translates. This N-point t-f space extends over the time 
interval from 0 to (N-l)T 2 and over the frequency interval from 
0 to (N-l) 3 in units of the normalized frequency fNT 4. 

5 

The iterated filter bank in FIG. 2 is used to generate the 
Wavelets which cover the t-f space in FIG. 1. Each filter stage 
5 consists of a high pass filter (HPF) and a low pass filter 
(LPF) . Output 6 of the LPF is subsampled by 2 which is 
10 equivalent to decimation by 2. This t-f space space is an N- 
dimensional complex vector metric space V 7. At stage m in the 
iterated filter bank, the remaining t-f space V m _ x 8 is 
partitioned into V m+1 9 and the Wavelet subspace W m+ i 10 . 

15 Scaling functions and Wavelets at each stage of this filter 

bank satisfy the following equations 



<p(n)= 2- ,/2 Zh q P(2n-q) V q 

(5) 

Kn)= 2- ,/2 Ig q ?>(2n-q) Vq 
q 

20 

where (p is the scaling function, \}/ is the Wavelet, HPF P 

coefficients are { Vq } , LPF P coefficients are { g q9 Vq } , and 

the equations apply to the stages 0,l,...,m-l. Identifying the 
scale parameter and using the previous Wavelet formulations 
25 enable these equations to be rewritten for stages p=0, 1, ...,m-l 
as 

^ (6) 
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For our application the HPF P and LPF P are quadrature mirror 
filters (QMF) with perfect reconstruction. This means they cover 
the subspace V p with flat responses over the subband frequency 
including the edges of the frequency subband, and the HPF P 
5 coefficients are the frequency translated coefficients for the 
LPF P : {g q =(-l) q h q , Vq} . 

Wavelet design using iterated filter bank starts with the 
selection of the scaling functions. Starting with a primative 
10 scaling function such as the one proposed by Daubechies ^one can 
use the iterated filter construction given by equations (5) and 
(6) to derive successive approximations to a desired scaling 
function which has properties that have been designed into it by 
the selection of the filter coefficients {g^Vq} at each level of 

15 iteration. The Wavelets can be derived from these scaling 
functions using the iterated filter construction or scaling 
equations (5) and (6) . 

Another use of the iterated filter construction is to 
20 design the scaling functions as Wavelets thereupon ending up with 
a larger set of Wavelets for multi-resolution analysis and 
synthesis as illustrated by Coifman's Wavelets in "Wavelet 
analysis and signal processing". 



SUMMARY OF THE INVENTION 

30 This invention for the design of multi-resolution Wavelet 

waveforms improves their performance for engineering and 
scientific applications, and in particular for applications to 
communications and radar. Current practice is to 1) design 
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traditional multi-resolution waveforms in the time domain using 
metrics which specify frequency performance and without 
consideration of Wavelet properties or to 2) design Wavelet 
multi-resolution waveforms using the time-domain iterated multi- 
5 resolution filtering approach with a set of scaling functions 
used to perform the filtering and without direct considerations 
of the frequency performance. These two separate design 

approaches yield fundamentally different waveforms. This 
invention provides a means to combine these two approaches to 
10 generate a new multi-resolution Wavelet waveform with the best 
properties of the traditional multi-resolution waveforms and the 
Wavelet multi-resolution waveform. 

This invention introduces innovations for the design of our 
15 Wavelet waveform that 1) include the ISI (intersymbol 
interference) , ACI (adjacent channel interference) , and QMF 
(quadrature mirror filter) requirements into the design algorithm 
along with non-linear modifications to the traditional passband 
and stopband frequency requirements which are used in the Remez- 
20 Exchange and eigenvalue LS algorithms, 2) select the frequency 
harmonics as the design coordinates, 3) use an appropriate 
subset of the available Fourier domain frequency harmonics as the 
design coordinates with the property that this subset of 
coordinates is a basis for the multi-resolution Wavelet waveform 
25 design, and 4) design the multi-resolution Wavelet waveform for 
no excess bandwidth ot=0. 

These innovations support the development of our new 
waveforms 1) that are generalizations of Wavelets in the 
30 frequency domain using a means which makes them useful for 
communications and radar over their t-f space, 2) that can be 
designed to be an orthonormal basis or orthonormal set of 
coordinates over the available frequency bandwidth with the 
implicit property that the excess bandwidth vanishes ct=0, to 



9 



within the accuracies allowed by communication and radar design 
implementations, and 3) that have multi-scale properties which 
allow a single dc waveform design to be used to uniquely define 
the complete set of waveforms over the t-f space for multi-scale 
5 applications Specific design algorithm examples developed in this 
invention disclosure are the LS design algorithms that use 
eigenvalue and gradient search techniques in the frequency domain 
to find the best waveform design which minimizes the 
corresponding LS error residual cost function that is a weighted 

10 linear sum of the residual error metrics. 

Wavelet waveform performance calculated in this invension 
disclosure from the application of these LS design algorithms, 
demonstrates the capability of this invention to provide a means 
to design waveforms that are improved over current practice. 

15 There are other algorithmic design algorithms which can be 
realized by modifications to this invention. An example given 
in this invention disclosure is the modification of these LS 
algorithms for application to constant amplitude bandwidth 
efficient (BEM) waveforms for communications. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS AND 
25 THE PERFORMANCE DATA 

The above-mentioned and other features, objects, design 
algorithms, and performance advantages of the present invention 
will become more apparent from the detailed description set forth 
30 below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 
elements, and in which: 
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FIG. 1 is a Wavelet N-point t-f space extending over the 
time interval (0, (N-l)T] and the frequency interval (0 , (N- 
1)/NT], which is tiled or covered by a set of orthonormal 
Wavelets at the scales p=0, 1, ...,m-l . 

5 

FIG. 2 is a Wavelet iterated filter bank used to generate 
the set of Wavelets which tile or cover the t-f space in FIG. 1. 

FIG. 3 illustrates the power spectral density (PSD) of our 
10 Wavelet waveform at dc, and the set of stopband and passband 
design requirements . 

FIG. 4 is the flow diagram of the LS metrics and cost 
15 functions and the final cost function used to find the optimal LS 
solution for our Wavelet waveform. 

FIG. 5 is flow diagram for the LS recursive solution 
algorithm to find the optimal harmonic coordinates that define 
20 the LS solution for our Wavelet waveform. 

FIG. 6 plots the dc PSD in dB units for the frequency 
response of our mother Wavelet at dc and the square-root raised- 
25 cosine (sq-rt r-c) waveforms with excess bandwidth parameter a= 
0.22, 0.40 vs. the normalized frequency in units of symbol rate 
1/T S . 

FIG. 7 plots the dc PSD in dB units for our waveforms vs. 
30 the normalized frequency fT s for our Wavelet waveform designed 
with modifications to the LS algorithms for application to 
constant amplitude BEM communications, and for a Gaussian 
minimum shift keying (GMSK) waveform. 
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FIG. 8 plots the amplitude of the dc radar ambiguity 
function vs. the normalized frequency fT p and normalized time tT c 
where the pulse time T p and chip time T c are both equal to the 
5 communications symbol interval T s for this example, for the new 
waveform and an unweighted chirp waveform. 
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DISCLOSURE OF THE INVENTION 

New Wavelet waveforms in this invention disclosure are 
generalizations of Wavelets in t-f space which enable them to be 

15 useful for communications and radar applications. This 
generalization is accomplished by 1) the introduction of a 
frequency translation, by 2) changing the orthonormality 
condition in equation (4) to apply to waveforms within the same 
space {q} and over the scales {p} with the inclusion of the 

20 frequency translation, and 3) by their characterization and 
design in the Fourier domain. With frequency translation the 
analytical formulation of these new waveforms as a function of 
the baseband or mother waveform centered at dc ( dc refers to the 
origin f=0 of the frequency space) becomes 

25 

V M M = 2" P V(2- p n -qM)e i2 ^ (p ^ T (7) 

30 

where f c (p,r) is the center frequency of the frequency translated 
dc waveform, at scale w p" and frequency index "r". The purpose 
of the frequency index "r" is to identify the center frequencies 
of the waveforms at the scale "p" in the t-f space. 
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These waveforms satisfy the complex orthonormality 
equations 



vi/* =1 iff p = kandq = m and r = v 
~ Y p,q,r Y k 9 m,v 

5 =0 otherwise (8) 

where is conjugation, and are generalizations of the 

orthonormality equations for the analytical Wavelets in (4) . 

10 The new waveforms in equation (7) expand the Wavelet 

analytical formulation to include a frequency variable. 
Wavelets are functions of the scale and translation parameters 
M p,q". The frequency variable together with the Fourier domain 
design are entirely new means for deriving these new waveforms 

15 as generalization of the traditional Wavelets 

These waveforms are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar in the t-f space. The basis vectors for 

20 this metric space V consist of a subset of the admissible set 
of scaled and translated waveforms { Pyq y<> Vp,q,r } derived from the 
dc waveform if/ as per equation (7) . An admissible waveform is 
any combination that covers F=t-f space. We are interested in 
the Fourier domain representation of the dc waveform (// in V , 

25 and in particular in a subset of the discrete Fourier transform 
(DFT) harmonic coefficients over the Fourier domain which we 
intend to use as the design coordinates. Starting with the z- 
transform and continuous Fourier transform, the DFT harmonic 
coefficients are defined by the following equations 

30 
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DFT Harmonic Coefficients 



y/(z) = ^y(n) z"" z-transform 

n 

y/(a)) = ^^(n) e' lfijn Fourier transform 

n 

y/ k = X^(n) DFT harmonic coefficients for V£ 

n 

k n 

= sin( (a)/2- riJNyN*)/ Nsin^ -akAT) 

k 

= ^y/ k [Harmonic interpolation for " k ?! ] 

k 

{k} = DFT frequency or harmonic coefficients such that f k = k/N f T 

where f k is the harmonic frequency corresponding to "k" 
N* = length of ^(n) 

where ^(co) = |^(6>)j 2 is the power spectral density of the dc 

waveform. These equations define the frequency representation of 
our new waveforms in terms of the available set of harmonic 
10 coefficients, which set is considerable larger than required for 
most applications . 

Our new waveforms are an orthonormal basis with no excess 
bandwidth, which properties are asymptoticly approached by our 

15 new waveforms to within design accuracies inherent in 
communications and radar. To demonstrate these properties we 
need to identify the structure of the dc waveform in V . We 
start with definitions for the parameters and coordinates in the 
following equation (10) . The waveforms derived from the dc 

20 waveform will be designed to be orthogonal over both time 
translates "MT" and frequency translates "1/MT" which 
respectively correspond to the Wavelet symbol spacing T S =MT and 
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the adjacent channel spacing 1/T S . This means the orthogonal 
spacing of the waveforms in V are at the time-frequency 
increments (MT, 1/MT) = (T s . 1/T S ) . In the interests of 
constructing our orthonormal multi-resolution waveforms to cover 
V it will be convenient to assume that M,L are powers of 2. We 
need the following definitions for the parameters and 
coordinates . 

Parameters and Coordinates (10) 

N' - Length of y which is an even function about the center and 
which spans an odd number of points or samples 
- ML + 1 where M, L are assumed to be even functions for 

convenience of this analysis 
= Number of points of \|/ 
M - Sampling interval for y/ 

= Spacing of y/ for orthogonality 
L - Length of y in units of the sample interval M 

= Stretching of v|/ over L sample intervals 
n = n 0 + 

= partitioning into an index n 0 over the sample 

length n 0 = 0,1 V ..,M- land an index n, = 0,1...,L-1 
over the sample intervals 
k =k 0 + k,L 

= partitioning into an index k 0 over the harmonic 
frequencies k 0 = 0,1,...,L-1 corresponding to the 
stretching and an index k, = 0,1,. ..,M over the harmonics 
frequencies corresponding to the admissible frequency 
slots for y/ 

The harmonic design coordinates are selected using the 
following observation. For most applications and in the following 
development it is assumed that the waveforms are spectrally 
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contained in the frequency interval 1/MT corresponding to the 
frequency spacing. This suggests the harmonic design coordinates 
be restricted to the subset of L harmonics { k 0 =0, 1, . . . , L-l } 
covering this spacing. These L harmonics correspond to the 
5 stretching of the mother waveform over the L repetition 
intervals . 

Obviously, for some applications as will be demonstrated 
later, the spectral containment is spread out over several 1/LT 
10 frequency increments whereupon one must increase the subset of 
design harmonics to possibly 2L, 3L or larger. 

, TheDFT equations for the dc waveform in (9) when 
rewritten in terms of the L harmonic design coordinates 

15 { ^ , Vk 0 } become: 

DFT equations for dc waveform 
20 if/ ko = ^ ^ (n) W^° n harmonic design coordinates 

n 

y/(n) =(1/N f ) £*F k W*° n new waveform defined in terms (11) 
of the L harmonic design coordinates {^ o ,Vk 0 } 

25 It will now be shown that the use of these L harmonic 

design coordinates is sufficient under time translates to be a 
basis for the corresponding subspace of V which means these 
waveforms provide a complete set of coordinates to describe this 
subspace. Wethe will use the theorems of Karhunen-Loeve and 

30 Mercer and will limit the demonstration to the dc waveform for 
simplicity and without loss of generality. We start by 
considering the expansion of a random complex sequence {z(n), 
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Vn} in a series of waveform coordinates consisting of time 
translates of the mother waveform \f/ . The sequence {z(n), Vn} 
is a zero-mean stationary random process which is orthonormal 
over the sample interval "M" and has a frequency spectrum which 
is flat and extends over the frequency range 1/LT which is 
centered at baseband corresponding to a zero frequency. This 
means the {z(n), Vn} cover the subspace of V corresponding to 
the scale of our dc waveform y/ and its time translates 

{ \f/(n-qM) = y/ q (ri), Vq } . In addition, the V is now considered to 

be extended over a time interval which is relatively large 
compared to the N-dimensional t-f space in FIG. 1,2 to avoid end- 
effects on the analysis. We start by approximating the sequence 
{z(n), Vn} by the {£(«), Vn} where: 



z(n)=ZZ q Kn-qM) 
q 

= ZZ q i^n 0 + (n,-q)M) 
q 

where the complex coefficien ts {z q } are 

derived from the original sequence using the (12) 
orthogonality properties of y/ 
Z q = Iz(n) ^*(n-qM) 

n 

= Iz(n)^(n) 

n 

The following equations prove that the coefficients {Z q , Vq} are 
orthonormal : 
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Z q Z q .* = Iz(An-qM)z*(An-q , MV q (An)^ q . *(An) (13) 

An 

= ^ Z|^q( n )| 2 since the sequence {z} is orthonormal 

n 

for qM time translates 
= 8 ' , with normalization of the energy of y/ 



5 Equations (12) and (13) together prove the Karhunen- 

Loeve's theorem which proves the following equation for the 
accuracy in approximating the stochastic sequence {z(n), Vn} by 
{z(w),Vn}. This accuracy is expressed by the expected "E(o)" 
squared error vx (o)" in this approximation: 

10 

E{z(n)-z(n)\ 2 }= 1 - J>(*o + M0f <"> 

We need to prove that the right hand side of this equation is 
15 zero which then proves that the approximating sequence is equal 
to the original sequence in the mean-square sense. In turn this 

proves that the new waveform coordinates { y/ g , Vq } are a basis for 
the original sequence {z(n), Vn} which is our goal. 

20 The right hand side of equation (14) when set equal to zero 

expresses Mercer's theorem so our goal is to prove Mercer's 
theorem. To do this we use the DFT of y/ in equation (11) and the 
coordinates in (10) to evaluate the right hand side of equation 
(14) . We find 

25 
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where T = (VL)^W*<* 

= sin(;r Ak 0 )/Zsin(;r Ak 0 1 L) 

= 1 for Ak 0 =0 (15) 
= 0 otherwise 
This proves that 

1. £k(" +n,M)\ 2 =0 Vn 0 



which proves that equation (14) reduces to 

5 

E[z(n)-z(n)\ 2 \= 0 (16) 

10 which as per the above proves that the set of multi-resolution 
waveforms is a basis. This proof easily generalizes to the 
multi-resolution waveforms at all of the scales {p} and time 
translates {q}, and to the expansion of the harmonic design 
coordinates over 2L, 3L, ... as required by the application. 

15 

Next we demonstrate that these waveforms have multi-scale 
properties, that we will demonstrate. We start with the 
observation that the Fourier domain design for the new waveforms 
provides a natural and easy way to derive the complete set of 

20 waveforms Vp,q} for the space V from the design of the dc 

waveform y/ by using the same invariant set of Fourier domain 
harmonic design coordinates {*P^,Vk 0 } derived for the dc multi- 
resolution waveform. This demonstration requires that we show 1) 
how the multi-scale transformations are implemented with the 
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design in the Fourier domain, and 2) how the waveform design 
remains invariant under scale changes. 

First consider the multi-scale transformation which derives 
5 the waveforms at the scale and shift parameters "p,q" from the dc 
waveform at scale "p=0" and centered at the origin "q=0". We 
begin by extending the parameters and coordinates in equation 
(10) , to include both scaling and subsampling or decimation in a 
form that is equivalent to the iterated filter bank construction 
10 which is used to derive the current Wavelet waveform using the 
filter scaling functions. Starting with the coordinates at scale 
xx p=0" the parameters and coordinates at scale "p=p" are given by 
the equations: 

15 

p=0 parameters and coordinates 

n = n 0 + ni M (17) 
n 0 - a 0 + ax 2 + ... + a m -i 2 ra " 1 
= M points 
20 M = 2 m 

ni = b 0 + bi 2 + ... +bi-i 2 1 " 1 

= L points spaced at M sample intervals 



25 p=p parameters and coordinates 

2' p n(l2 p ) = scaled by "2" p " 

and subsampled or decimated by 2 P :1 

= n 0 (p) + ni(p)M 

2' p n 0 (i2 p ) = scaled by "2" p " 
30 and subsampled or decimated by 2 P :1 

= n 0 (p) 

= a p + a p _i2 + ... + a p+m _i 2 p+m 1 

= M points spaced at 2 P sample intervals 
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2" p n 1 (il p ) = scaled by "2' p " 

and subsampled or decimated by 2 P :1 
= n x (p) 

= b p + b p+1 2 + ... + bp^.x 2 p+1 - x 
5 = L points spaced at M2 P sample intervals 

together with the observation that the sampling interval W T" is 
increased to "2 P T" under the scale change from "p=0" to "p=p" and 
subsampling or decimation from "1:1" to "2 P :1". Combining these 
10 equations with the analytical formulation in (7) and the Fourier 
domain representation in (11) enables the waveforms at the 
parameters "p,q" to be written as a function of the Fourier 
domain harmonic design coordinates: 
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¥ M M =2- (V2 K2- p n-qM)e^ (p * T 

= (2- p/1 /AT) 2>, o w^ (n(p)HqM) e 



iln f c {ps)n{p)2*T (18) 



for all admissible scale, translation, and frequency index 
20 parameters "p,q,r". 

Next we need to demonstrate that the frequency domain 
design in (11) remains invariant for all parameter changes and in 
particular for all scale changes. This multi-scale property 

25 expresses the accordion behavior of the design in that the 
Wavelets at different scales are simply the stretched and 
compressed versions of the mother waveform with the appropriate 
frequency translation indices. This multi-scale invariancy means 
that the design for a M=16 channel filter bank remains the same 

30 for M=100 or M=10,000 channel filter banks, when the overlap L 
and the performance goals remain constant. To demonstrate this 
invariant property across scales, we consider the multi- 
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resolution waveform at scale "p" with the other parameters set 
equal to zero for convenience "q=0, r=0" and without loss of 
generality. The Fourier domain frequency response can be 

evaluated starting with the original formulation in equation (7) : 



DFT at "p, q=0, r=0' 



(19) 



k 0 n(p) 



-(JN'2 p T -k 0 )n(p) 



sm(7T (fN'2 p T-k 0 )) 
N'sin(;r (fN2 p T-k 0 )/N') 



= [Harmonic interpolation for n k 0 "] 
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This only differs from the harmonic representation in equation 
(9) in the restriction of the design coordinates to the subset of 

harmonic coefficients { y/ k ,Vk 0 } and the stretching of the time 

15 interval to "2 P T" corresponding to the scale "p=p" . The harmonic 
interpolation functions are observed to remain invariant over 
scale changes upon observing that the frequency scales as 
"f~l/2 p T" which means the "f requency*time" product remains 
invariant with scale changes as per the fundamental property of 

20 the waveforms. This means the harmonic interpolation functions 
remain invariant with scale change and therefore the frequency 
response remains an invariant. This demonstrates the waveform 
design is an invariant across the waveform scales which means we 
only need a single design for all scales or resolutions of 

25 interest. 
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LS design algorithms for new waveform will be described to 
illustrate the advantages our waveform has over current designs. 
The two LS algorithms described are the eigenvalue and the 
gradient search which respectively can be reduced to algorithms 
5 which are equivalent to current eigenvalue and Remez-exchange 
waveform design algorithms for application to a uniform filter 
bank. We consider the t-f space which is spanned by a uniform 
polyphase filter bank consisting of M channels at the frequency 
spacing f w =l/MT where T is the digital sampling interval, and 

10 the filter waveform FIR time response is stretched over L 
sampling time intervals T s . This polyphase filter bank is ideally 
decimated which means the filter ouput sample rate 1/T S is equal 
to the channel-to-channel spacing T S =MT, equivalent to stating 
that there is no excess bandwidth a=0 . Our design for this 

15 topology is immediately applicable to an arbitrary set of multi- 
resolution filters through the scaling equation (18) which gives 
the design of our waveform at arbitrary scales in terms of our 
design of the dc waveform. 

20 For this polyphase filter bank used to construct the dc 

waveform or filter impulse response, our LS example design 
algorithms will use 5 metrics consisting of the 2 prior art 
passband and stopand metrics, and the 3 new metrics consisting 
of the ISI, ACI, and QMF, and solve the LS minimization problem 

25 using as design coordinates the subset of harmonic coordinates 
which are a basis. Since our two example LS design algorithms 
only differ in the use of an eigenvalue LS optimization and the 
use of a gradient search LS optimization, the flow diagrams for 
the construction of the cost functions and the solution for the 

30 optimal waveform will be identical. However, there are 
differences in the construction of the cost functions from the 
respective metrics and in the iterative solution mathematics. 
Both LS solutions for the harmonic design coordinates minimize 
the weighted sum of the error residuals or cost functions from 
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the 5 metrics. These design coordinates are the Fourier harmonics 

{y/ k ,Vk 0 } for {k 0 = 0, 1, L-l} . Resulting algorithms are easily 

extended to the applications requiring the design coordinates to 
cover 2L, 3L, ... harmonics. 

Frequency domain design coordinates are related to the 
waveform time domain digital samples or coordinates as follows. 

Mappings of frequency to time (20) 

Time domain design coordinates {^(«), Vn}are real and 
symmetric and can be represented by the reduced set {h t (n), 
n=0, l / ... / ML/2} 

h t (n) = ^(0) for n=0 

= 2y/{n) for n=l,2,...,ML/2 

= time domain design coordinates 

Frequency domain harmonic design coordinates {y/ k ,Vk 0 } are 

real and symmetric and can be represented by the reduced 
set {h f (k), k=0,l,,...,L-l} 

h f (k) = y/ k ^ for k=k 0 =0 

= 2y/ ko for k=k 0 =l, 2, L-l 

= frequency domain design coordinates 

Mapping of the frequency coordinates {h f (k), k=0, 1,...,L-1} 
into the time coordinates {h t (n), n=0, 1, ML/2 } is defined 
by the matrix transformation to within a scale factor 
h t = B h f 

where 

h f = (hf(0), ...,h f (ML/2) ) c transpose of column vector 
h t = (h t (0), ... f h t (ML/2) ) t transpose of column vector 
B = (ML/2 +1) x L matrix 
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= [ B kn ] matrix of row k and column n elements B kn 
B kn = 1 for n=l 

= 2 cos(27ckn/ ML) otherwise 



5 

wherein the N'=ML+1 has been replaced by ML since a single end 
point has been added to the the FIR to make it symmetrical for 
ease of implementation for the example Wavelet being considered 
10 with a sample at the mid-point that makes the number of samples 
N' an odd number. 

Passband and stopband metrics and cost functions are 
derived with the aid of FIG. 3 which defines the po PSD 

15 parameters of interest for the passband and stopband of the PSD 
*P(<y) for communications applications. Requirements for radar 
applications include these listed for communications. Referring 
to FIG. 3 the passband 11 of the waveform PSD is centered at dc 
(f=0) since we are designing the dc or baseband waveform, and 

20 extends over the frequency range (Dp extending from -(Op/2 to +a>p/2 
12 in units of the radian frequency variable co=27ifT 13. The 
frequency space extends over the range of f=-l/2T to f=+l/2T 
which is the frequency range in FIG 1 translated by -1/2T so 
that the dc waveform is at the center of the frequency band. 

25 Quality of the PSD over the passband is expressed by the passband 
ripple 14. Stopband 15 starts at the edge 16 of the 
passbands of the adjacent channels +/-oo a /2 16 and extends to 
the edge of the frequency band a>=+/-7i 17 respectively. 
Stopband attenuation 18 at +/-co a /2 measures the PSD isolation 

30 between the edge of the passband for the dc waveform and the 
start of the passband for the adjacent channels centered at +/-go s 
19. Rolloff 20 of the stopband is required to mitigate the 
spillover of the channels other than the adjacent channels, onto 
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the dc channel. Deadband or transition band 21 is the interval 
between the passbands of contiguous channels, and is illustrated 
in FIG. 3 by the interval from Op/2 to G) a /2 between the dc 
channel and adjacent channel at a> a . Waveform sample rate co s 22 
5 is the waveform repetition rate. For the LS example algorithms, 
the waveform sample rate is equal to the channel-to-channel 
spacing for zero excess bandwidth. Therefore, 1/T S = (d s /2jiT=1/MT 
which can be solved to give (o 3 =2n/M for the radian frequency 
sampling rate of the filter bank which is identical to the 
10 waveform repetition rate. 

We start by rewriting the DFT equations for the dc waveform 
in (11) as a function of the {h f (k), k=0, 1, L-l} 

y/{o)) - ^^(n)cos(nty) 

n 

15 = c'Bh f using (21) and the definition of the vectorV 1 (21) 

c = ( 1 , cos(tf)), cos( (ML/2) 6)) y transpose of column vector 



which is equivalent to the equation for ys(a)) in (9) expressed as 
a linear function of the {h f (k), k=0, 1, L-l } . An ideal "c" 
20 vector "c r " will be introduced for the passband and the stopband 
in FIG. 3 in order to identify the error residual 8y/((D) at the 
frequency "<s>" in meeting the ideal passband and stopband 
requirements. The ideal PSD is flat and equal to "1" for the 
passband, and equal to "0" for the stopband. We find 

25 

Error residuals for passband and stopband (22) 

c r = (1, 1,..., l) t passband ideal "c" 
30 = (0,0,. ..,())* stopband ideal "c" 

Sc = c r - c error vector 
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= Residual error in meeting the ideal spectrum at 11 o>" 



The LS metric for the passband and stopband can now be 
constructed as follows for the eigenvalue and the LS optimization 
(or equivalehtly , the LS algorithm) design algorithms 



Passband and stopband metrics 



(23) 



J (band) = * J|£^(tf>)| 2 do) Eigenvalue 



band 



Eigenvalue 
LS 



where 



band = [0, C0p ) 

= (<0 S , 7t] 
1 

R = 

band 

5i|/ 



IMI - 



passband 
stopband 

— |(B t && t B)d^ 

band 

= L x L matrix 
= ( 8v|/((o 1 ) 5h/(g) u ) ) t 
= vector of error residuals at the 

frequencies ©i , o>u across the band 

norm or length of the vector (o) and which 

includes a cost function for the errors of 
the individual components 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 
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coordinates {h f (k), Jc=0,l,...,L-l} whereas with the LS approach it 
is sufficient to give the LS metrics as vector norms with 
imbedded cost functions or an equivalent formulation. 

5 QMF metrics express the requirements on the deadband that 

the PSD's from the contiguous channels in FIG. 3 add to unity 
across the deadband [a>p , co s ] in order that the filters be QMF 
filters. By suitable modification of the error vector 8c, the 
previous construction of the passband and stopband metrics can be 
10 modified to apply to the deadband. We find 



Deadband metrics (24) 



15 J (deadband) = ||<^(<y)| do Eigenvalue 

deadband dead ^ and 

(24) continued 
= h' f R h f Eigenvalue 

= ||<V || 2 LS 

where 

20 6c = c r - c(a>) - c(n/M - co) 



where c (co) =c as defined in (22) and (23), and c(7t/M - a>) = c at 
the offset frequency "7t/M - co" corresponding to the overlap of 
25 the contiguous filters over the deadband. 

Orthonormality metrics measure how close we are able to 
designing the set of waveforms to be orthonormal over the t-f 
space, with the closeness given by the ISI and the ACI . ISI and 
30 ACI errors are fundamentally caused by different mechanisms and 
therefore have separate metrics and weights to specify their 
relative importance to the overall sum of pfre L$ metrics. ISI 
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is a measure of the non-orthogonality between the stream of 
waveforms within a channel as per the construction in FIG. 3. 
On the other hand, ACI is a measure of the non-orthogonality 
between the waveform within a channel and the other waveforms in 
5 adjacent channels. This means the stopband performance metric 
has a significant impact on the ACI due to the sharp rolloff in 
frequency of the adjacent channel, and the ACI metric is then a 
measure of the residual non-orthogonality due to the inability of 
the stopband rolloff in frequency from completely eliminating the 
10 ACI errors. 

We assume that the received waveform is identical to the 
filter waveforms and is transmitted at the filter output sample 
intervals equal to MT seconds. The second assumption means we 
15 are assuming the receiver is synchronized with the received 
signal. Since there is no information lost by sampling 
asynchronously with the received waveform, we are free to make 
this synchronization assumption without loss of generality. ISI 
metrics are derived in the following set of equations. 

20 

ISI metrics (25) 

Mapping of h f into \\f 
25 y = (-ML/2) , ... , V|/ (ML/2)) t transpose of column vector 

= H h f 

H = (ML+1) x L matrix of elements H kn 

Hkn = 1 for n=l 

= 0.5 cos (27ikn/ML+l) otherwise 
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Offset matrix A 

A = L x (ML+1) matrix of elements Akn 

Akn = 0 for row vector k=l 
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= [0 0 ... vj/(-ML/2) ... V|/ ( (L-k)M+l) ] row vector k 
12 kM LM+1 column elements 



ISI error vector 8E 

8E = L x 1 column vector 



= A H h 



f 



ISI metric 

J (ISI) = 8E fc 8E Eigenvalue 
10 = Non-linear quadratic function of h f 

= Ue\\ 2 ls 



ACI metrics are derived using the ISI metric equations with the 
15 following modifications . 

ACI metrics (26) 

Mapping of h f into vy is the same as developed for ISI 
20 Offset matrix A elements are changed as follows to apply to 

Channel 1: 

A kn = 0 for row vector k=l 

= [0 0 ... hM-ML/2)^ ... v|/( (L-k)M+l) W { J'~ k) ] 
1 2 kM+1 LM+1 

25 which means the ACI error vector 8E is 

8E = L x 1 column vector 
= A H h f 

ACI metric for the two contiguous channels 
30 J (ISI) = 2 SE* 1 8E Eigenvalue 

= Non-linear quadratic function of h f 

= 2 \\SE\\ 2 LS 
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where the factor "2" takes into account there are two contiguous 
channels or one on either side of the reference channel 0 in FIG. 
3. Because of the fast rolloff of the frequency spectrum the 
5 addition of more channels into the ACI metric is not considered 
necessary, although the functional form of the ACI metric in (26) 
allows an obvious extension to any number of adjacent channels 
which could contribute to the ACI - 

10 Cost function J for the LS algorithms is the weighted sum 

of the LS metrics derived in (23) , (24) , (25) , (26) . The LS 
algorithms minimize J by selecting the optimal set of frequency 
coordinates {h f (k),Vk} for the selected set of parameters used 
to specify the characteristics of the dc waveform, frequency 

15 design coordinates, LS metrics, and weights. Cost function and 
optimization techniques are given by the equations 



Cost function J 



(27) 
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J 




weighted sum of the LS metrics J (metrics) 



where 



metrics = passband, stopband, deadband, ISI, ACI 
{w (metrics) , ^/metrics} = set of weights 



25 




= 1 normalization 



Optimization goal 

Goal: minimize J with respect to the selection of 



the {h f (k),Vk} 
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Optimization algorithms 

Two algorithms are the Eigenvalue and the LS optimization 



where the eigenvalue optimization algorithm uses the non- 
linear quadratic formulations of the LS metrics and the LS 
optimization algorithm uses the norm formulations for the 
LS metrics. 



FIG. 4 is a summary of the LS metrics and the construction 
of the cost function J. Design parameters 23 are L, N f =ML+l, 
number of DFT samples per FIR length n_fft for implementation of 

10 the LS algorithms, passband radian frequency (Dp, stopband radiam 
frequency <o a , waveform repetition rate in radian frequency <d s , 
selection of the set of design coordinates {h f } to be used in 
the optimization, and the metric weights {w (metrics )} . Output 
parameters are the set of harmonic design coordinates {h f } that 

15 minimize J. Band metrics 24 are the passband, stopband, and 
deadband metrics defined in equations (23) , (23) , (24) 
respectively. Interference metrics 25 are the ISI and ACI 
metrics defined in equations (25) and (26) respectively. LS cost 
function J 26 is the weighted linear sum of the metrics defined 

20 for the band 24 and the interference 25 as defined in 
equation (27) . 

FIG. 5 is a flow diagram of the LS recursive solution 
algorithm. There are two loops with the topology constructed so 

25 that the outer loop 27 is an iteration over the set of metric 
weights {w (metrics) } , and the inner nested loop 28 is the 
recursive or iteratived LS solution to find the optimal {h f } for 
the given design parameters and weights where each step refers to 
an iteration step in the solution. A recursive LS solution is 

30 required to be able to solve the highly non-linear interference 
metrics as well as the band metrics when one chooses to use non- 
linear techniques to construct these metrics. 
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Consider the inner nested loop 28. The recursive LS 
solution starts with the initial step i=0 29 which begins with 
the selection of the input design parameters and weights and the 
selection of an initial set of values for the {h f } 30, Next the 
5 band metrics are calculated 31. Since this is the initial step 
i=0 32 the cost function J is restricted to the linear band 
metrics in order to fine the linear approximation to {h f } 33 to 
initialize the non-linear solution starting with step i=l 34 
wherein the highly non-linear interference metrics are included 

10 in J. Following the signal flow, for step i=l the band metrics 
31 and interference metrics 35 are calculated, weighted, and 
summed to form the cost function J 36. An iterative solution 
algorithm 37 finds the best approximation in step i=l to the 
{h f } which minimizes J using the approximation to {h f } from the 

15 previous step i=0 to linearize the search algorithm coefficients 
for the eigenvalue and LS step i=l iteration. If there is no 
convergence to the correct solution for {h f } this recursive 
sequence of calculations is repeated for step i=2 34. This 
recursive solution technique is repeated for subsequent steps 

20 until there is convergence 38 whereupon one exits this inner 
loop . 

Consider the outer loop. After exiting the inner loop 38 
the solution for {h f } is tested to see if it meets the 

25 performance goals 39. If not, a new set of metric weights is 
selected and the inner loop is initialized i=0 40 and the 
inner loop is used to find the next solution for {h f } . This 
process continues until the performance goals are met or 
adequated approximated 39 whereupon the algorithm is exited 

30 with the final solution set {h f } 41. 
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10 



FIG. 4,5 derive the new multi-resolution Wavelet \y at 
baseband. Scaling equations (7) , (18) can be used to derive the 
corresponding digital Wavelet at scales p,q,k with parameters 
M,L wherein the digital frequency parameter "k" replaces the 
frequency index "r". Equation (7) can be rewritten to define the 
scaling equation which derives this Wavelet V|/ P/q , M ,L,k from v|/: 



New Wavelet scaled from mother Wavelet (28) 

M>p, q ,M,L,k = M> P ,q,k with parameters M,L 

= 2" p/2 vj/(2" p n-qM) exp(j2rc k 2" p n /ML) 



15 

There are at least two methods of scaling used to implement 
equation (28) . In method 1 corresponding to the standard Wavelet 
assumption, the subsampling or decimation of n by the factor 2 A p 
stretches the time between contiguous samples by the factor 2 A p 

20 while keeping constant the number of samples M per Wavelet \\f 
sample interval, and therefore reduces the frequency band by the 
factor 2 A p. For a uniform filter bank, this means the filter 
frequency bands are reduced by the factor 2 A p. In method 2, the 
sampling is held constant corresponding to a constant frequency 

25 band assumption. The factor 2 A (-p) is divided out in the Wavelet 
expression which is equivalent to the M being scaled to 2 A pM 
corresponding to increasing the number of samples per Wavelet 
interval by the factor 2 A p. For a uniform filter bank, this 
means method 2 corresponds to increasing the maximum number of 

30 channels by the factor 2 A p over the same frequency band. 



Applications of this new invention to both communications 
and radar will be given using these example algorithms and other 
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algorithms supported by this invention. These new waveforms are 
considered for the applications: l)to replace the square-root 
raised cosine waveform (sq-root rc) which is extensively used for 
the third generation (3G) CDMA communications, 2) to replace the 
5 Gaussian minimum shift keying (GMSK) waveform for constant 
amplitude bandwidth efficient (BEM) applications, and 3) as a 
candidate waveform for synthetic aperture radar (SAR) and real 
aperture radar (RAR) applications. 

10 CDMA communications application for the current and the new 

3G CDMA considers a waveform designed with this new invention as 
a possible replacement for the sq-rt rc waveform with bandwidth 
expansion parameter a = 0.22 to a = 0.40 . This notation means that 
for a = 0.22 the spectral efficiency is (symbol rate/bandwidth) = 

15 1/1+a =1/1.22 = 0.82 = 82%. A basic advantage of the waveform 
is the potential for a symbol rate increase within the same 
bandwidth with an increase in the spectral efficiency to «100% 
depending on the application and operational constraints. The dc 
power spectral density or power spectrum (PSD) of the waveform is 

20 compared to the PSD for the sq-rt r-c in FIG. 6. Plotted are the 
measured PSD in dB 42 versus the frequency offset from dc 
expressed in units of the symbol rate 43 Plotted against the 
normalized frequency offset are the dc PSD for the new waveform 
44, the sq-rt r-c with a=0.22 45, and the sq-rt r-c with 

25 a=0.40 46. It is observed that the PSD for the new waveform 
rolls off faster that that for the sq-rt r-c which means that the 
new mr waveform will support an increased symbol rate for a given 
available frequency band while satisfying the inherent 
requirements for low ISI and MAI (multiple access interference) . 

30 

Constant amplitude BEM application of the new waveform 
indicates that it is a viable candidate for replacing the current 
preferred modulation waveform which is the GMSK. The GMSK finds 
applications for transmitters which operate their HPA(s) 
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amplifiers in a saturation mode in order to maximize their 
radiated power from the HPA(s), and which require a BEM PSD to 
avoid excessive spreading of the transmitted power. Simulation 
data for the dc PSD is plotted in FIG. 7 for the new waveform BEM 
5 and the GMSK. Plotted are the measured dc PSD in dB 47 versus 
the frequency offset from dc expressed in units of the bit rate 
48. Plotted against the normalized frequency offset are the dc 
PSD for the new waveform BEM 49 and the GMSK 50 , for a length 
parameter L=10 where L is the length of the phase pulse in terms 
10 of the phase pulse repetition rate. The significance of this 
example data is that the new waveform has the potential to be 
designed to offer a PSD which is less spread out than the current 
GMSK, and therefore an improved BEM waveform. 

15 Radar RAR and SAR application of the waveform indicates 

that it is a viable candidate to replace the current chirp 
waveforms for wideband signal transmission, when combined with 
pseudo-random phase codes. Results of the simulation for the new 
waveform and an unweighted frequency chirp waveform are given in 

20 FIG. 8. Plotted are the ambiguity function for the new waveform 
51 and the unweighted frequency chirp waveform 52. The dc 2- 
dimensional radar ambiguity function 53 is plotted as a function 
of the frequency offset in units of fT p and the time offset in 
units of t/T c where T p is the phase-coded radar pulse length or 

25 length of the phase code and T c is the phase code chip length. 
The chip length is identical to the waveform repetition interval 
T 3 so that T c = T s . It is observed that the new waveform has the 
potential for significant improvements in the ambiguity function 
and by implication in the performance. 

30 

Preferred embodiments in the previous description are 
provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
35 and the generic principles defined herein may be applied to other 
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embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 
embodiments shown herein andis to be accorded the wider scope 
consistent with the principles and novel features disclosed 
herein. 
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5 ABSTRACT OF THE DISCLOSURE 

ABSTRACT 

The prcocnt invention dcocribco a now method for the design 
e£ — multi - resolution — waveformo — fee — improve — system — performance, 
dcoign — flexibility, applicability, aad — reduce — the — cooto — e# 

10 implementation, compared — to — the — current — ene^r-. £ke — multi - 
rcoolution waveform designs arc a generalization of Wavelets to 
the — Fourier domain. — Thio — generalization uses — dcoign algorithmo 
which include application metrico to improve oyotcm performance, 
allowo a oinglc dcoign to bo uocd for all ocales, and allowo 

15 each deoign to be characterized by a relatively few parameters. 

Preliminary oimulationo indicate that — waveform — aftd filter 

performance — can be — significantly — improved — compared — fee — current 
techniques — — communications — applications — te — code — division 
multiple — acccoo — (CDMA) — waveformo — and — filters, — communicationo 

20 conotant amplitude — bandwidth — efficient — waveformo, aftd — radar 

oyotcm applicationo . 

Wav elets for multi resolution waveformo and filters ar e 
de signed to meet specific application requirements, arc defin ed 

25 by a relatively small number of frequency harmonics, arc com^e^ 
irft cludc a frequency tranolation, and can be designed to -be 
orthogonal with no exceso bandwidth for communica -^jrefte 
ap plicationo. Frequency dcoign harmonico and frequency 
^Tr anslation capability enable the waveformo at multip jre 

30 resolutions to be derived from a single design by scaling— 

de lation, translation, and frequency tranolation parameter s 
Be oign is illustrated by Matlab 5.0 code to generate a linear 
#*- ltcr waveform using an iterative least squares eigcnval rae 
ap proach to minimize the non linear least squareo coot functi en-r 



ee d to scale thio waveform dcoign for multi resolution 
ap plication opacified by tho dilation, time translation, — a**dt 
£g equcncy translation parameters. Additional results arc giv e** 
-ge- r a constant amplitude minimum shift - keying bandwidth - efficient 
5 modulation waveform and for a synthetic aperture radar waveform. 

A method for designing Wavelets for communications and 
radar which combines requirements for Wavelets and finite impulse 
response FIR filters including no excess bandwidth, linear 

10 performance metrics for passband, stopband, quadrature mirror 
filter QMF properties, intersymbol interference, and adjacent 
channel interference, polystatic filter design requirements, and 
non-linear metrics for bandwidth efficient modulation BEM and 
synthetic aperture radar SAR. Demonstrated linear design 

15 methodology finds the best design coordinates to minimize the 
weighted sum of the contributing LS error metrics for the 
respective performance requirements. Design coordinates are 
mapped into the optimum FIR symbol time response. Harmonic 
design coordinates provide multi-resolution properties and enable 

20 a single design to generate Wavelets for arbitrary parameters 
which include dilation, down-sampling, up-sampling, time 
translation, frequency translation, sample rate, symbol rate, 
symbol length, and set of design harmonics. Non-linear 
applications introduce additional constraints. Performance 

2 5 examples are linear communications, BEM, and SAR. 
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radar which combines requirements for Wavelets and finite impulse 
response FIR filters including no excess bandwidth, linear 

10 performance metrics for passband, stopband, quadrature mirror 
filter QMF properties, intersymbol interference, and adjacent 
channel interference, polystatic filter design requirements, and 
non-linear metrics for bandwidth efficient modulation BEM and 
synthetic aperture radar SAR. Demonstrated linear design 

15 methodology finds the best design coordinates to minimize the 
weighted sum of the contributing LS error metrics for the 
respective performance requirements. Design coordinates are 
mapped into the optimum FIR symbol time response. Harmonic 
design coordinates provide multi-resolution properties and enable 

20 a single design to generate Wavelets for arbitrary parameters 
which include dilation, down-sampling, up-sampling, time 
translation, frequency translation, sample rate, symbol rate, 
symbol length, and set of design harmonics. Non-linear 
applications introduce additional constraints. Performance 

25 examples are linear communications, BEM, and SAR. 
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APPLICATION NO. 09/826,118 

TITLE OF INVENTION: Wavelet Multi-Resolution Waveforms 
INVENTOR: Urbain A. von der Embse 



Currently amended CLAIMS 
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APPLICATION NO. 09/829,118 

INVENTION: N|^j Multi-Resolution Waveforms 

INVENTORS: Urbain Alfred A . von der Embse 



CLAIMS 

WHAT IS CLAIMED IS: 

Claim 1. ( currently amcndcd deleted ) A means for the dcoign 
of new multi - resolution wavcformo and filters in the Fourier 
domain with properties which 

provide — extensions — of the Wavelet — concept to — fcke — Fourier 
domain or cquivalontly the frequency domain 

provide single waveform — designs — for all — of the waveform s eet 

multiple scales 

provide designs which requires a relatively few design 

coordinates compared to the FIR time response samples for the 
wavcf orm 

provide design methodologies — which can incorporate application 

metrics to improve the waveform performance 
provide designs that allow the use of direct design methodologies 
that circumvent the need to solve a Wavelet iterated filter bank 
construction to obtain the waveforms thereby providing improved 
flexibility to meet the application goals 

Claim 2. (deleted) ^; A — means — §€H? — the — design — — multi - 
resolution — wave forms — aftd — filters — ift — fe-he — frequency — domain — en? 
cquivalcntly the Fourier of the t - f space for waveform and filter 
applications, — with properties which 

provide a means for single design to be used for all of the 

waveforms at — the multiple resolutions or frequency bands 
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provide a means for designs that require a relatively few 
design — coordinates — compared — fce — fcke — number — e4 — digital 
oamplco covered by the Wavelet 

provide a means for the designs to include frequency aftd 

5 time application metrics to improve the waveform performance 

provide a means to use direct design methodologies in the 
frequency - time domain c for all of the waveforms and filters for 
the multiple frequency bands 

10 

Claim 3. ( currently amondcd deleted ) 

A method for the design of multi - resolution waveforms which 

allows — Fourier domain techniques — to be used. Properties 

can include — some or more of the listed properties in — (4-)- 
15 and (2) 



Claim 4. (deleted) 



20 A method for the design of multi resolution waveforms which 

eeft — incorporate — Fourier — domain — techniques — into — design 
methodologies — which — eaf* — include — analytical — aftd — iterated 
filter bank construction design techniques. 



25 

Claim 5. ( currently amcndcd deleted ) — 

A method for the analysis and design of multi resolution 
waveforms using Fourier domain techniques which talec 
advantage of the new invention disclosures on the 
30 characterization and design of multi-resolution waveforms 

in the — Fourier domain. 
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Claim 6. ( currently amcndod deleted ) A f*ew 

formulation — if* — multi - resolution — waveform — ae — a 

function of the dc multi - resolution waveform which addo tho 
concept — — a — frequency — index — that — allows — the — multi - 
5 rcoolution waveform to be placed arbitrarily throughout the 

£ — space — thereby — — avoiding — the — restrictions — e4 — the 
Wavelet — iterated — filter — construction — §ef — tiling — a — 
space, and 2) — allowing the new multi rcsoloution waveforms 

tre — he — used — for multi resolution — communications aftd £e*r 

10- bandwidth on - demand communications application, — in place of 

traditional Wavelets. 



15 

Claim 7. (new) An iterative eigenvalue least-squares LS 
method for designing digital mother Wavelets at baseband for 
means for the design of new multi-resolution waveforms and 
filters, said method comprising steps: 
20 power spectral density PSD representative requirements for said 
mother Wavelet \\f frequency (0 response \\f & {a>) in a multi- 
channel filter bank, specify 

a) passband frequency range for waveform transmission, 

b) stopband spacing between adjacent filters, 
25 c) bounds on ripple over said passband, 

d) stopband filter attenuation, 

e) rolloff with frequency outside stopband, 

f) quadrature mirror filters QMF require the sum of said 
PSD' s for contiguous filter responses to be flat over 

30 deadband which is said stopband, 

q) symbol-to-symbol interference ISI, 

h) adjacent channel interference ACI, 
said LS error metrics to measure said requirements (a) -(h) are 
derived as functions of said Wavelet i|/(n) assuming 
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i) T and 1/T are sample interval and sample rate equal to 

Nyquist sample rate, 

j) \\f is real and symmetric about n=0, 

k) n= 0, + /-1, . . +/-ML/2 digital index over said \y, 

5 1) M is interval between contiguous said \\f , 

m) 1/MT is said \\f symbol rate and channel-to-channel 

separation, 

n) L is length of said \\f in units of said M, 

said multiple-resolution properties require said LS metrics 

10 to be constructed as functions of said Wavelet Fourier 

harmonics \\f k (k) with k=Q, +/-!,. - . , + /-(N k -l) and 

N k >L is a design parameter, 

it is sufficient to use positive n=0,l,. . ML/2 and k=0,l,. . 

. ,N k -l since said \y (n) and \\f k (k) are real and symmetric, 

15 ML/2+1 x N k matrix bw wherein "x" reads "by" maps i|/ k (k) into 



v|/(n) to within a scale factor by equations 



vi/ (n) = 


Zk bw(n+l,k+l)\j7 k (k) for 


n>0 , k>0 , 


bw (n+1, 


k+1) = 1 


for n=0, 




= 2 cos (2n nk/ML) 


otherwise, 



20 = row n+1, column k+1 element of bw, 

said LS error metrics are converted by said bw mapping into 

quadratic forms in said h k equal to J(band)=h k 'R h k wherein 
said h k ' is the transpose of h k and said R is a real square 
symmetric matrix of LS errors in meeting said requirements, 
25 LS ISI,ACI error metrics J (ISI ) , J (AC I ) are derived as non-linear 
quadratic forms in h and converted by said bw matrix to the 
non-linear quadratic form in h k equal to J ( I S I ) =5E ' 8E , 
J(ACI) =28E'5E wherein SE =AHh k is a column vector and 
matrix "A" in the matrix product AH is a function of said h 
30 hereby introducing said non-linearity, and said AH differ 

for ISI and ACI error metrics, 
LS cost function J is the weighted sum of said LS error metrics 
J = Z w(LS metric) J(LS metric) 
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with summation over said LS metrics= passband, stopband, 
QMF deadband, ISI, ACI with normalized weights 
Zw(LS metric) =1, 
said weights are free design parameters, 
5 said iterative eigenvalue LS algorithm at each step finds the 

optimum eigenvalue and eigenvector which minimize said 

quadratic form J in h k for a constant said "A", 

said eigenvector is the optimum h k which minimizes said J and 

said bw equation derives the corresponding optimum h which 

10 minimizes said J, 

step 1 in said iterative algorithm finds said optimum eigenvalue, 
eigenvector, h k , h of J reduced by deleting said non-linear 
ISI and ACI LS quadratic error metrics, 
said h is used to evaluate said "A" matrices for step 2, 
step 2 finds said optimum eigenvalue, eigenvector, h k , h for 

minimum J using said "A" from step 1, 

said h is used to evaluate said "A" for step 3, 
steps 3,4, etc. continue until said minimum J converges to a 
steady value and, 

said optimum i|/ k (k) uses said bw to calculate optimum \j/(n) for 

implementation as said Wavelet FIR digital waveform and 

filter time response. 



15 



20 



25 Claim 8. (new) An LS method for designing digital mother 

Wavelets at baseband for multi-resolution waveforms and filters, 
said method comprising steps: 

said PSD waveform representative requirements and assumptions 

are recited in (a)-(n) in Claim 7, 

30 said multiple-resolution properties require said LS metrics 

to be constructed as functions of said v|/ k (k), 

said LS error metrics for said passband, stopband, and QMF 

deadband requirements are derived as squared vector norm 

functions of said h and converted by said bw matrix into 

5 



J (band) = 



- wherein l|Bh^| is the vector norm of the 



column vector Bh k and said B is the matrix of LS errors in 



meeting said requirements and wherein said squared vector 

norm is suitable for LS optimization, 
5 LS ISI,ACI error metrics J(ISI) , J(ACI) are derived as squared 

vector norm functions equal to J(ISI= ||5e|| 2 , J(ACI)=2 ||5e|| 2 

using said column vectors 5E =AHh k in claim 7 , 
LS cost function J is said weighted sum of said LS error metrics 

equal to J=Zw(LS metric) J (LS metric) defined in claim 7, 

10 an LS gradient search algorithm finds optimum h k (k) to 
minimize J, 

step 1 of said LS gradient search algorithm uses a Remez- 

exchange algorithm to find said optimum h k (k) for said J 

reduced to said passband and stopband LS metrics, 

15 step 2 uses the estimated h k (k) from step 1 to initialize said 
gradient search, 

step 3 selects one of several available gradient search 

algorithms, gradient search parameters, and stopping rules, 

step 4 implements said , algorithm, parameters, and stopping rule 
20 selected in step 3 to derive said optimum h k (k) to minimize 

J and, 

said optimum h k uses said bw to calculate optimum \j/(n) for 

implementation as said Wavelet FIR digital waveform and 
filter time response. 

25 

Claim 9. (new) Wherein said mother Wavelet generates 
multi-resolution dilated Wavelets, comprising steps and design: 
said Wavelet parameters are 
30 a) scaling parameter p dilates sampling by factor 2 A p 

equivalent to sub-sampling by factor 2^p, 



6 



b) translation parameter q translates said \y by qM digital 
samples, 

c) frequency offset k is set by design, 

d) symbol repetition interval said M remains constant, 
5 e) Wavelet length said L in units of said M remains 

constant, 

step 1 uses said design harmonics \\f k to generate said FIR time 

response \j/(n p ) at baseband with said bw equation 

V(n p ) = Z k bw(n p +l,k+l) \|/ k (k) 

10 recited in claim 7 with 

bw = 2 cos (2n n p k/ML) for n p >0 
wherein n p =n/2 A p is n sub-sampled or equivalently dilated 
by the factor 2 A p, 

step 2 uses said i|Mn p ) to construct said multi-resolution 

15 Wavelet vjJ p ,o, M .L.k with equation 

*&>,q,M,L,k = 2 A (-p/2) vj/(n p -qM) exp(j2rc kn p /ML) 

which is said FIR time response for parameters p,q,M,L,k 

wherein the subset p,M,L are the scale parameters, 

design of said multi-resolution Wavelet includes 

20 f) said T for n is increased to T2 A p for n p / 

g) said 1/T is reduced to l/T2 A p, 

h) said y symbol rate 1/MT equal to said channel-to-channel 
separation is reduced to l/MT2 A p in Hz and, 

i) said \\f length (ML+1)T in seconds is stretched to 
25 (ML+l)T2 A p in seconds. 



Claim 10 (new) Wherein said mother Wavelet generates multi- 
resolution constant sample rate dilated Wavelets, comprising 
30 steps and design: 

said Wavelet parameters are 

a) said p dilates said \y to increase said length from 
ML+1 to M P L+1 where M p =M2 A p is the dilated interval 
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b) 


said 


q 


translates said \\? by qM p digital samples, 


c) 


said 


k 


is set by design, 


d) 


said 


M p 


=M2 A p is dilated M, 


e) 


said 


L 


remains constant, 



step 1 uses said design harmonics i[/ k to generate said FIR time 

response v|/(n p ) at baseband with said bw equation 

V(n) = I k bw(n+l,k+l) v|/ k (k) 

recited in claim 7 with 

10 bw = 2 cos (27i nk/MpL) for n>0, 

step 2 uses said Vj/(n) to construct said multi-resolution Wavelet 
U£p,q,M,L,k with equation 

v^p, q ,M,L,k = 2 A (-p/2) H/(n-qM p ) exp(j27t kn/M p L) 

which is said FIR time response for parameters p,q,M,L,k, 

15 design of said multi-resolution Wavelet includes 

f) said T remains constant, 

g) said 1/T remains constant, 

h) said \\f symbol rate 1/MT equal to said channel-to-channel 
separation is reduced to 1/M P T=1/MT2 A p in Hz and, 

20 i) said \\f length (ML+1)T in seconds is stretched to 
(ML2 A p+l)T in seconds. 



Claim 11. (new) Wherein said mother Wavelet generates multi- 
25 resolution up-sampled Wavelets, comprising steps and design: 
said Wavelet parameters are 

a) said p up-samples said digital sampling rate 1/T to 2 A p/T, 

b) said q translates said vy by qM digital samples, 
c) said k is a design parameter, 

30 d) said M is constant, 

e) said L is constant, 
step 1 uses said design harmonics h f ,v|; k to generate said FIR 

time response v|/(n p ) at baseband with said equation 
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V|Mn p ) = Z k bw(n p +l, k+1) \y k (k) 



recited in claim 7 with 

bw = 2 cos (27i ripk/ML) for n p >0 
wherein n p is n up-sampled by the factor 2 A p and defined 
5 by equations 

n p = n p + n 2 A p 

n p = 0, 1,2, . . . , 2 A p-l 

wherein n p is the index over the additional samples added 

to each sample n by said up-sampling, 

10 step 2 uses said \|/(n p ) to construct said multi-resolution Wavelet 
\(/p,q,M,L,k with equation 

y^P, q ,M,L,k = 2 A (-p/2) v(n p -qM) exp(j2ft kn p /ML) 

which is said FIR time response for parameters p,q,M,L,k, 

design of said multi-resolution Wavelet includes 
15 f) said T is decreased to T/2 A p, 

g) said 1/T is increased to 2 A p/T, 

h) said \\f symbol rate 1/MT equal to said channel-to-channel 
separation is increased to 2 A p/MT in Hz and, 

i) said \|/ length (ML+1)T in seconds is reduced to 
20 (ML +l)T/2 A p in seconds. 



Claim 12 (new) Wherein said multi-resolution Wavelets have 

25 properties comprising : 

said scale parameters p,M,L and said design parameter 1/T specify 

said multi- resolution Wavelets at baseband and said q, k 

specify time, frequency translations from baseband, 

said design harmonics \|/ k (k) of mother Wavelet are said design 

30 coordinates for multi-resolution Wavelets, 

said design harmonics y k {k) use said bw matrix to generate said 
multi-resolution Wavelet baseband time response vy(n) for said 
dilation, dilation of Wavelet length, and up-sampling as 
recited in Claims 9-11 and which is translated 
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in time and frequency to said multi-resolution Wavelet 

said design harmonics are few in number compared to said 

Vj/(n) , 

5 said \y is designed to support a bandwidth-time product B f T=l+g 
with no zero excess bandwidth g=0, 

said multi-resolution Wavelets are designed to behave like an 

accordion in that at different scales said Wavelets are 

stretched and compressed versions of the mother Wavelet with 

10 appropriate time and frequency translation, 

said optimization techniques in claims 7,8 assume said i|/(n) 

symmetric about n=0 and are applicable to other arrangements 
of vj/(n) with self-evident modifications, 
optimization algorithms for finding said optimum set of i|/ k (k) 
15 use said linear LS waveform and filter design methods 

recited in claims 7,8 and also use other methods and, 
said linear waveform and filter LS design methods can be modified 

to design waveforms for applications including bandwidth 

efficient modulation BEM and synthetic aperture 
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APPLICATION NO. 09/826,118 

TITLE OF INVENTION: Wavelet Multi-Resolution Waveforms 
INVENTOR: Urbain A. von der Embse 



Clean version of how the CLAIMS will read. 
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APPLICATION NO. 09/829,118 

INVENTION: Multi-Resolution Waveforms 

INVENTORS: Urbain A. von der Embse 



CLAIMS 
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WHAT IS CLAIMED IS: 



Claim 1. (deleted) 

Claim 2. (deleted) 

Claim 4. (deleted) 

Claim 5. (deleted_ 

15 Claim 6. (deleted) 

Claim 7. (new) An iterative eigenvalue least-squares LS 
method for designing digital mother Wavelets at baseband for 
multi-resolution waveforms and filters, said method comprising 
20 steps: 

power spectral density PSD representative requirements for said 
mother Wavelet y frequency a> response i|/©((0) in a multi- 
channel filter bank, specify 

a) passband frequency range for waveform transmission, 
25 b) stopband spacing between adjacent filters, 

c) bounds on ripple over said passband, 

d) stopband filter attenuation, 

e) rolloff with frequency outside stopband, 

f) quadrature mirror filters QMF require the sum of said 
30 PSD's for contiguous filter responses to be flat over 

deadband which is said stopband, 

g) symbol-to-symbol interference ISI, 

h) adjacent channel interference ACI, 

said LS error metrics to measure said requirements (a) -(h) are 
35 derived as functions of said Wavelet i|/(n) assuming 
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i) T and 1/T are sample interval and sample rate equal to 

Nyquist sample rate, 
j) v(/ is real and symmetric about n=0, 

k) n= 0, . . +/-ML/2 digital index over said 

5 1) M is interval between contiguous said \\f, 

m) 1/MT is said \\f symbol rate and channel-to-channel 
separation, 

n) L is length of said V|/ in units of said M, 
said multiple-resolution properties require said LS metrics 
10 to be constructed as functions of said Wavelet Fourier 

harmonics vj/ k ( k) with k=0, . . +/-(N k -l) and 

N k >L is a design parameter, 
it is sufficient to use positive n=0,l,. . . , ML/2 and k=0,l,. . 
. ,N k -l since said v|/(n) and y k (k) are real and symmetric, 
15 ML/2+1 x N k matrix bw wherein "x" reads xx by" maps \|/ k (k) into 
i|/(n) to within a scale factor by equations 
v|/(n) = Zk bw(n+l,k+l)\j/ k (k) for n>0, k>0, 
bw(n+l,k+l) = 1 for n=0, 

= 2 cos(2tc nk/ML) otherwise, 
20 = row n+1, column k+1 element of bw, 

said LS error metrics are converted by said bw mapping into 

quadratic forms in said h k equal to J(band)=h k 'R h k wherein 
said h k ' is the transpose of h k and said R is a real square 
symmetric matrix of LS errors in meeting said requirements, 
25 LS ISI,ACI error metrics J(ISI) , J(ACI) are derived as non-linear 
quadratic forms in h and converted by said bw matrix to the 
non-linear quadratic form in h k equal to J(ISI) =5E'8E, 
J(ACI)=25E'5E wherein 8E =AHh k is a column vector and 
matrix "A" in the matrix product AH is a function of said h 
30 hereby introducing said non-linearity, and said AH differ 

for ISI and ACI error metrics, 
LS cost function J is the weighted sum of said LS error metrics 
J = Z w(LS metric) J(LS metric) 
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with summation over said LS metrics= passband, stopband, 
QMF deadband, ISI, ACI with normalized weights 
Zw(LS metric) =1, 
said weights are free design parameters, 
5 said iterative eigenvalue LS algorithm at each step finds the 
optimum eigenvalue and eigenvector which minimize said 
quadratic form J in h k for a constant said "A", 
said eigenvector is the optimum h k which minimizes said J and 

said bw equation derives the corresponding optimum h which 
10 minimizes said J, 

step 1 in said iterative algorithm finds said optimum eigenvalue, 
eigenvector, h k/ h of J reduced by deleting said non-linear 
ISI and ACI LS quadratic error metrics, 
said h is used to evaluate said "A" matrices for step 2, 
15 step 2 finds said optimum eigenvalue, eigenvector, h k/ h for 
minimum J using said "A" from step 1, 
said h is used to evaluate said "A" for step 3, 
steps 3,4, etc. continue until said minimum J converges to a 
steady value and, 
20 said optimum \\t k (k) uses said bw to calculate optimum y(n) for 
implementation as said Wavelet FIR digital waveform and 
filter time response. 



25 Claim 8. (new) An LS method for designing digital mother 

Wavelets at baseband for multi-resolution waveforms and filters, 
said method comprising steps: 

said PSD waveform representative requirements and assumptions 
are recited in (a)-(n) in Claim 7, 
30 said multiple-resolution properties require said LS metrics 
to be constructed as functions of said V|/ k (k), 
said LS error metrics for said passband, stopband, and QMF 

deadband requirements are derived as squared vector norm 
functions of said h and converted by said bw matrix into 
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J (band) = Bh 

IE 



2 wherein Bh jJ is the vector norm of the 



column vector Bh k and said B is the matrix of LS errors in 
meeting said requirements and wherein said squared vector 
norm is suitable for LS optimization, 
5 LS ISI,ACI error metrics J ( ISI) , J (ACI ) are derived as squared 

vector norm functions equal to J (ISI=||5e|| 2 , J (ACI ) =2 ||5e|| 2 

using said column vectors 8E =AHh k in claim 7, 
LS cost function J is said weighted sum of said LS error metrics 
equal to J=Xw(LS metric) J (LS metric) defined in claim 7, 
10 an LS gradient search algorithm finds optimum h k (k) to 
minimize J, 

step 1 of said LS gradient search algorithm uses a Remez- 

exchange algorithm to find said optimum h k (k) for said J 
reduced to said passband and stopband LS metrics, 
15 step 2 uses the estimated h k (k) from step 1 to initialize said 
gradient search, 
step 3 selects one of several available gradient search 

algorithms, gradient search parameters, and stopping rules, 
step 4 implements said algorithm, parameters, and stopping rule 
20 selected in step 3 to derive said optimum h k (k) to minimize 

J and, 

said optimum h k uses said bw to calculate optimum \|/(n) for 

implementation as said Wavelet FIR digital waveform and 
filter time response. 

25 



Claim 9. (new) Wherein said mother Wavelet generates multi- 
resolution dilated Wavelets, comprising steps and design: 
said Wavelet parameters are 
30 a) scaling parameter p dilates sampling by factor 2^p 

equivalent to sub-sampling by factor 2 A p, 
b) translation parameter q translates said y by qM digital 
samples, 
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c) frequency offset k is set by design, 

d) symbol repetition interval said M remains constant, 

e) Wavelet length said L in units of said M remains 
constant, 

step 1 uses said design harmonics v|/ k to generate said FIR time 
response v|/(n p ) at baseband with said bw equation 
i|/(n p ) = Z k bw(n p +l,k+l) V|/ k (k) 
recited in claim 7 with 

bw = 2 cos (27i n p k/ML) for n p >0 
wherein n p =n/2 A p is n sub-sampled or equivalently dilated 
by the factor 2 A p, 
step 2 uses said i)/(n p ) to construct said multi-resolution 
Wavelet \\f p , q ,M f L,k with equation 

M> P ,q,M,L,k= 2 A (-p/2) v|/(n p -qM) exp(j2rc kn p /ML) 
which is said FIR time response for parameters p,q,M,L,k 
wherein the subset p,M,L are the scale parameters, 
design of said multi-resolution Wavelet includes 

f) said T for n is increased to T2 A p for n p , 

g) said 1/T is reduced to l/T2 A p, 

h) said \|/ symbol rate 1/MT equal to said channel-to-channel 
separation is reduced to l/MT2 A p in Hz and, 

i) said i|/ length (ML+1)T in seconds is stretched to 
(ML+1) T2 A p in seconds . 

Claim 10 (new) Wherein said mother Wavelet generates multi- 
resolution constant sample rate dilated Wavelets, comprising 
steps and design: 
said Wavelet parameters are 

a) said p dilates said \y to increase said length from 
ML+1 to M P L+1 where M p =M2 A p is the dilated interval 
between contiguous i|/'s, 

b) said q translates said V|/ by qM p digital samples, 



c) said k is set by design, 

d) said M p =M2 A p is dilated M, 

e) said L remains constant, 

step 1 uses said design harmonics \|/ k to generate said FIR time 
response V|/(n p ) at baseband with said bw equation 
v|Mn) - Zk bw(n+l,k+l) v|/ k (k) 
recited in claim 7 with 

bw = 2 cos (2ti nk/MpL) for n>0, 
step 2 uses said vj/(n) to construct said multi-resolution Wavelet 
H/ P ,q,M,L,k with equation 

V P , q ,M,L,k= 2 A (-p/2) \|/(n-qM p ) exp(j2* kn/M p L) 
which is said FIR time response for parameters p,q,M,L,k, 
design of said multi-resolution Wavelet includes 

f) said T remains constant, 

g) said 1/T remains constant, 

h) said \j/ symbol rate 1/MT equal to said channel-to-channel 
separation is reduced to l/M p T=l/MT2 A p in Hz and, 

i) said \\f length (ML+1)T in seconds is stretched to 
(ML2"p+l)T in seconds. 

Claim 11. (new) Wherein said mother Wavelet generates 
multi-resolution up-sampled Wavelets, comprising steps and 
design: 

said Wavelet parameters are 

a) said p up-samples said digital sampling rate 1/T to 2 A p/T, 

b) said q translates said \|/ by qM digital samples, 

c) said k is a design parameter, 

d) said M is constant, 

e) said L is constant, 

step 1 uses said design harmonics h f ,i|/ k to generate said FIR 
time response vj/(n p ) at baseband with said equation 
v|/(n p ) = Zk bw(n p +l,k+l) \|/ k (k) 



recited in claim 7 with 

bw = 2 cos(27c ripk/ML) for n p >0 
wherein n p is n up-sampled by the factor 2 A p and defined 
by equations 

n p = n_p + n 2 A p 
n_p =0,1,2,. . . , 2 A p-l 
wherein n_p is the index over the additional samples added 
to each sample n by said up-sampling, 
step 2 uses said \y (n p ) to construct said multi-resolution Wavelet 
y P , q ,M,L,k with equation 

Vp, q ,M,L f k= 2 A (-p/2) v|/(n p -qM) exp(j2n kn p /ML) 
which is said FIR time response for parameters p,q,M,L,k, 
design of said multi-resolution Wavelet includes 

f) said T is decreased to T/2 A p, 

g) said 1/T is increased to 2 A p/T, 

h) said \\f symbol rate 1/MT equal to said channel-to-channel 
separation is increased to 2 A p/MT in Hz and, 

i) said \\f length (ML+1)T in seconds is reduced to 
(ML +l)T/2 A p in seconds. 



Claim 12 (new) Wherein said multi-resolution Wavelets have 
properties comprising : 

said scale parameters p,M,L and said design parameter 1/T specify 
said multi- resolution Wavelets at baseband and said q,k 
specify time, frequency translations from baseband, 

said design harmonics V|/ k (k) of mother Wavelet are said design 
coordinates for multi-resolution Wavelets, 

said design harmonics \}/ k (k) use said bw matrix to generate said 
multi-resolution Wavelet baseband time response V|/(n) for 
said dilation, dilation of Wavelet length, and up- 
sampling as recited in Claims 9-11 and which is translated 



in time and frequency to said multi-resolution Wavelet 

Vp,q,M,L,k/ 

said design harmonics vj/ k (k) are few in number compared to said 
vj/(n) , 

5 said \\j is designed to support a bandwidth-time product B f T=l+a 
with no zero excess bandwidth a=0, 
said multi-resolution Wavelets are designed to behave like an 
accordion in that at different scales said Wavelets are 
stretched and compressed versions of the mother Wavelet with 
10 appropriate time and frequency translation, 

said optimization techniques in claims 7,8 assume said vj/(n) 

symmetric about n=0 and are applicable to other arrangements 
of vy(n) with self-evident modifications, 
optimization algorithms for finding said optimum set of V|/ k (k) 
15 use said linear LS waveform and filter design methods 

recited in claims 7,8 and also use other methods and, 
said linear waveform and filter LS design methods can be modified 
to design waveforms for applications including bandwidth 
efficient modulation BEM and synthetic aperture 
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